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A theoretical overview of the phenomenon of spontaneous magnon decays in quantum antifer- 
romagnets is presented. The intrinsic zero-temperature damping of magnons in quantum spin 
systems is a fascinating many-body effect, which has recently attracted significant attention in 
view of its possible observation in neutron-scattering experiments. An introduction to the the- 
ory of magnon interactions and a discussion of necessary symmetry and kinematic conditions for 
spontaneous decays are provided. Various parallels with the decays of anharmonic phonons and 
excitations in superfluid 4 He are extensively used. Three principal cases of spontaneous magnon 
decays are considered: field-induced decays in Heisenberg antiferromagnets, zero-field decays in 
spiral antiferromagnets, and triplon decays in quantum-disordered magnets. Analytical results are 
compared with available numerical data and prospective materials for experimental observation 
of the decay-related effects are briefly discussed. 
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I. INTRODUCTION 

According to conventional wisdom a quasiparticle is 
presumed well-defined until proven not to be. The text- 
book picture of magnons as long-lived excitations weakly 
interacting with each other works very well for many 
magnetic materials. Nevertheless, in a number of recently 
studied spin systems magnons neither interact weakly nor 
remain well-defined even at zero temperature. This is due 
to the so-called spontaneous quasiparticle decay, a spec- 
tacular quantum- mechanical many-body effect, which is 
the topic of the current Colloquium. 

Magnetic excitations have been the subject of inten- 
sive theoretical and experimental studies since the first 
half of the last century. Historically, the concept of spin 
waves wa s introduced by |Bloch| ( |193Q[ ) and, indepen- 
dently, by Slater (1930) for the wavelike propagation of 



spi n flips in ordered ferromagnets; see also Ho ddeson et 
at1(|1987|). Already in his earliest work|Bloch| (|193Q|) de- 



scribed spin waves as weakly interacting excitations obey- 
i ng Bo se statistics. Subsequently, Hols tein and Primakoff| 
(1940) introduced second quantization of spin waves in 



* Electronic address: 
t Electronic address: 



mike.zhitomirsky@cea.fr 



sasha@uci.edu 



terms of bosonic operators, the approach widely used to 
this day. The first use of the word "magnon," which con- 
cisely unifies the notions of an elementary quasiparticle 
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and a magnet ic quantum, can be traced back to |Pomer- 
anchuk (|1941|), who ascribe d this term to Lev Landau, 



Walker and Slack (1970). Although the terms "spin 



waves" and "magnons" are sometimes used to distinguish 
between long- and short- wavelength excitations, respec- 
tively, in the following we use both names interchange- 
ably. 

In addition, the fundamental contribution to the the- 
ory of spin waves in magnet ic ins ulators was made by | An- 



derson! jKuto^ and [Dyson| ( |1956a|bD . These 



and subsequent studies developed a comprehensive pic- 
ture of spin waves in conventional quantum ferromagnets 
and antiferromagnets, which can be found in various re- 



views 


and books, see Van Kranendonk and Van Vleck 


(1958 


), Akhiezer et al (|1968D, Mattis_dl981 ), White 


(1983 


), Borovik- Romanov and Sinha (1988), Manousakis 


(1991 


), and Majlis (2007). Invented in the 1950s, inelas- 



tic neutron-scattering spectroscopy ( |Brockhouse , 1995) 
has quickly become a major experimental method of in- 
vestigating spin waves and other magnetic excitations in 
solids. With a variety of specialized techniques currently 
used at numerous neutron facilities around the world, 
it remains the most powerful and versatile research tool 
for studies of magnetic spectra. Other microscopic tech- 
niques to investigate excitations in magnetic materials 
are electron spin resonance dGurevich and Melkov[ |l996 ) 
and light scatteri ng methods, such a s resonant inelastic 
X-ray scattering ( Beaurepaire , 2006). 



Usually, the magnon-magnon interactions play a minor 
role at low temperatures. For example, as T — >• devia- 
tion of the spontaneous magnetization of a ferromagnet 
from its saturated value follows Bloch's law AM ex T 3 / 2 , 
the result obtained by considering magnons as noninter- 
acting bosons. Interactions yield a correction to Bloch's 
law that scales as T 4 , which is still smaller than the "kine- 
matic" correction proportional to T 5 / 2 due to deviation 



of magnon dispersion from its parabolic form (Dyson 



1956b). In antiferromagnets, effects of magnon-magnon 
interaction survive at zero temperature due to zero-point 
motion in the ground state. An estimate of their role is 
given by the magnitude of the quantum correction to 
the magnon dispersion law obtained in the harmonic, 
or semiclassical, approximation. Even in low dimension 
(2D) and even in the "extreme quantum" case of S = 1/2, 
such a correction for the square-lattice Heisenberg anti- 
ferromagnet amounts to a mere y 16% rescaling of the 
harmonic dispersion (Manousakis, 1991). 



Recently, there has been a growing body of theoretical 
and experimental studies showing that the conventional 
picture of magnons as weakly interacting quasiparticles is 
not always correct. The first class of systems exhibiting 
significant deviations from the conventional semiclassical 
dynamics is weakly coupled antiferromagnetic chains, for 
which the spin-wave approximation is not a good start- 
ing point. Their high-energy spectra are well-described 
by propagating one-dimensional spi nons, the natural ba- 



sis for excitations in spin-^ chains (Coldea et al 



Kohno et al, 2007 Schulz, 1996). Since the development 
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of a quantitative theory for the ID to 3D crossover in 
weakly coupled chains is still in progress, we shall not 
discuss this physics here. 

The second class is, in a sense, simpler and more puz- 
zling at the same time. It is comprised of quantum mag- 
nets with well-ordered ground states whose dynamical 
properties are yet quite different from the predictions of 
harmonic spin- wave theory ( Zheng et al\ 2006a). In the 
following, we focus on this second class of unconventional 
magnetic systems. It includes quantum magnets with 
noncollinear spin ordering. The key role in their unusual 
dynamical properties is played by cubic anharmonicities, 
or three-magnon processes, which originate from the non- 
collinear ity of spins. Three-magnon processes are ab- 
sent in conventional collinear antiferromagnets, making 
their spin-wave excitations intrinsically weakly interact- 
ing. The most striking effect of three-magnon interac- 
tions is the finite lifetime of spin waves: not only at 
nonzero temperatures when magnon linewidth is deter- 
mined by scattering off thermal excitations, but also at 
T = when damping is produced by spontaneous magnon 
decay. 

In order to put the phenomenon of spontaneous 
magnon decays in a broader context, we note that two 
other subfields of condensed matter physics also deal with 
decays of bosonic excitations. These are anharmonic 
crystals (Kosevich , 2005 ) and superfluid liquids and gases 
QLifshitz and Pitaevsknf|1980| |Pethick and Smith] |2008[ ) . 



Cubic anharmonicities are naturally present in these sys- 
tems, but their role differs both quantitatively and qual- 
itatively from that in quantum spin systems. In crystals, 
the anharmonic corrections are generally small due to 
small mean fluctuations of ions around equilibrium posi- 
tions, the condition satisfied for the majority of solids be- 
low the Debye temperature. In superfluids, the strength 
of cubic interactions is determined by the condensate 
fraction and may again become small, as is the case of 
superfluid 4 He. Quantum magnets are unique in that (i) 
cubic anharmonicities may be tuned on and off by an 
external magnetic field or by changing the lattice geome- 
try, and (ii) the existing variety of magnetic compounds 
provides diverse forms of the magnon dispersion law 
which plays key role in spontaneous quasiparticle decay. 
Still, various parallels and analogies between the three 
branches of condensed matter physics prove to be helpful 
and will be extensively used in the subsequent discussion. 

The rest of the article is organized as follows. Section[TI| 
contains a general discussion of the origin of cubic as well 
as quartic anharmonicities in quantum magnets. The en- 
ergy and momentum conservation conditions for spon- 
taneous decays and the decay-induced singularities in 



two-dimensional spin systems are analyzed in Sec. Ill In 
Sec. |IV| we provide a description of spontaneous magnon 
decays in the square-lattice Heisenberg antiferromagnet 
in strong external field, while Sec. |V| deals with another 
prototypical noncollinear magnet: the triangular-lattice 
Heisenberg antiferromagnet in zero field. A few general 
scenarios of the decay of long- wavelength excitations are 
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considered in Sec. I VII Section rVTU is devoted to the dis- 
cussion of the role of cubic anharmonicities in gapped 
quantum spin liquids. Finally, a summary and outlook 
are presented in Sec. |VIII 



As pointed out by Anderson (1984), the symmetry- 



II. MAGNON-MAGNON INTERACTIONS 

In the classical picture, the interaction between 
magnons is related to the amplitude of spin precession 
and can become arbitrarily small as the amplitude de- 
creases. Zero-point motion of spins puts a lower bound 
on the amplitude and leads to non-negligible interaction 
effects in quantum magnets. Usually, the consideration 
of magnon-magnon interactions begins with a discussion 
of bosonic representations for spin operators. Here we in- 
tentionally skip this step, postponing technical details to 
the subsequent sections. Instead, we discuss the general 
structure of effective bosonic Hamiltonians defined solely 
by the symmetry of the magnetic ground state and cor- 
responding excitations. 



A. Four-magnon interactions 

To draw an analogy between different bosonic systems 
we begin with the Hamiltonian of a normal Bose gas 



Ho 



k 



(1) 



7 X) y 4 0) (ki , k 2 ; k 3 , k 4 ) a ^ a£ a « k3 a k 



k t 



where is the kinetic energy, V^°^(ki, k^; k3, k^) is the 
two-particle scattering amplitude, the ellipsis stands for 
n-particle interactions, and momentum conservation is 
assumed from now on for various k-summations. The in- 
teraction term in Eq. conserves the number of parti- 
cles, a natural constraint for liquid helium or cold atomic 
gases. At the same time, an effective bosonic Hamilto- 
nian with exactly the same structure as in Eq. ([!]) de- 
scribes the quantum Heisenberg ferromagnet (| Hoist em] 



and Primakoff , 1940 Oguchi I960). Although the num- 



ber of quasiparticles may not be conserved, in a ferro- 
magnet the conservation is enforced by the invariance 
of the ground state under an arbitrary rotation about 
the magnetization direction. Because of this symmetry 
the ground state and all excitations are characterized by 
definite values of the z-projection of the total spin S z , 
from S z ot = NS and down. Consequently, the matrix 
elements of the Hamiltonian vanish for states with dif- 
ferent spin projections, resulting in the particlelike form 
of the interaction term in Eq. ([!]). In other words, in 
the Heisenberg ferromagnet every magnon has an intrin- 
sic quantum number AS Z = —1, which is conserved in 
magnon scattering processes. 

The quasiparticle-number conservation does not hold 
for interaction processes in quantum antiferromagnets. 



broken antiferromagnetic ground state corresponds to a 
superposition of states with different values of the total 
spin. Hence, elementary excitations cannot be assigned 
with a definite value of spin: the spin of the spin wave 
ceases to exist. Because of that, the effective Hamilto- 
nian of a quantum antiferromagnet contains additional 
interaction terms, 

^4 = ^X F 4 (1) (ki,k 2 ,k 3 ;k4) (4 1 4 2 4 3 a k 4 + h - c -) 

ki 

1 ^y 4 (2) (k 1 ,k 2 ,k 3 ,k 4 ) (4 1 4 a 4 B 4 4 + h - c -)> ( 2 ) 



4! 



that do not conserve the 
for example, |Harris et al. 



number of excitations; see, 
first term in 



^ fll97ip . The 
Eq. ([2| describes decay and recombination processes of 
one magnon into three and vice versa, whereas the sec- 
ond, so-called source term, corresponds to creation and 
annihilation of four particles out of (into) an antiferro- 
magnetic vacuum. 



B. Three-magnon interactions 

Another type of anharmonicity that may appear in an 
effective bosonic Hamiltonian is the three-particle inter- 
action term. In a generic quantum system with noncon- 
served number of particles anharmonicities of all orders 
are present, beginning with the lowest-order cubic terms, 
which couple one- and two-particle states. For instance, 
cubic terms represent the dominant anharmonicity for 
lattice vibrations (Kosevich 2005). However, additional 



symmetry restrictions may forbid such cubic interactions 
in antiferromagnets. 

In the case of a collinear antiferromagnet, the remain- 
ing SO (2) rotational symmetry about the Neel vector 
direction (z axis) prohibits cubic terms. Although the 
single-magnon state is not an eigenstate of the S z ot op- 
erator, such a state still preserves an odd parity under 
the 7r rotation about the z axis. In contrast, the two- 
magnon states are invariant under this symmetry op- 
eration. Consequently, the coupling between one- and 
two-particle sectors is strictly forbidden in collinear anti- 
ferromagnets and cubic terms do not occur in their ef- 
fective bosonic description. In hindsight, this lack of 
low-degree anharmonicities explains why the spin-wave 
theory works quantitatively well even for spin-^ systems: 
magnets with collinear spin structures have excitations 
that are intrinsically weakly-interacting. 

On the other hand, three-magnon interaction terms 
must be present if the spin-rotational symmetry is com- 
pletely broken in the ground state, i.e., when the mag- 
netic structure becomes noncollinear. This is realized, 
for example, due to spin canting in an applied magnetic 
field or because of competing interactions in frustrated 
antiferromagnets. The general form of the three-boson 



interaction is given by 

^3 = ^E^3 (1) (ki,k 2 ;k 3 ) «at 2 a k3 + h.c.) 
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where the first term describes the two-particle decay 
and recombination processes while the second is another 
source term. 

There exists a deep analogy between the noncollinear 
quantum magnets and supernuid boson systems regard- 
ing the presence of cubic interactions in their effective 
Hamiltonians. The Bose gas in the normal state is in- 
variant under the group of gauge transformations U(l) 
and is described by the Hamiltonian ([!]). The macro- 
scopic occupation of the lowest-energy state (ak=o) 7^ 
below the supernuid transition leads to breaking of the 
U(l) symmetry and, at the same time, to the appearance 
of the cubic terms by virtue of the Bogolyubov substitu- 
tion 



4 1 4 2 a ks a k 4 4X2^3^0) + • ■ ■ ■ 



t J 



(4) 



In other words, the vacuum state with the Bose conden- 
sate absorbs or emits an extra particle, enabling noncon- 
serving processes between excitations. 

An analogous consideration of the noncollinear mag- 
netically ordered state is possible, leading to the same 
qualitative answer. In fact, it was realized a long time ago 
by |Matsubara and Matsuda] (|1956[) and by |Batyev and 



Braginskii (1984) that there is a one-to-one correspon- 
dence between breaking the SO (2) rotational invariance 
about the field direction for Heisenberg and planar mag- 
nets and breaking the U(l) gauge symmetry at the su- 
pernuid transition. Noncollinear magnetic structures sta- 
bilized by competing exchange interactions in zero field 
belong to a different class as they spontaneously break 
the full SO (3) rotational symmetry without the help of 
an external field. Still, the absence of any remaining 
symmetry constraints for their elementary excitations, 
apart from energy and momentum conservation, permits 
all possible anharmonic terms including the cubic ones. 

It is necessary to mention that cubic anharmonicities 
due to the long-range dipolar interactions have been dis- 
cussed already in the early works on the quantum the- 



ory of spin waves in ferromagnets (Akhiezer et al t 1968 
Holstein and Primakoff[ |194Q[ ) where they play an impor- 



tant role in the low-frequency ma gnetization dynamics 
flChernyshev] |2Q12| [Gurevich and Melkov| |1996[ ). How- 
ever, the effect of cubic terms is completely negligible at 
higher energies because of the smallness of the dipole- 
dipole interaction compared to the exchange energy. In 
contrast, cubic anharmonicities discussed in this work 
are of exchange origin and have a profound effect on the 
magnon excitation spectra. 

Anharmonic interactions are not restricted to cubic 
and quartic processes. The usual spin- wave expansion 
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FIG. 1 (color online) . Zero-temperature self-energy diagrams 
due to the (a) and (b) three-magnon, and (c) and (d) four- 
magnon interactions. 



based on the Holstein-Primakoff transformation for spin 
operators produces an infinite series of interaction terms. 
Typically, these higher-orders terms are small due to the 
spin-wave expansion parameter (a ( a)/2S <C 1 associated 
with each order of expansion. We disregard them in the 
following analysis as they do not alter our results, either 
qualitatively or quantitatively. 



C. Magnon self-energies 

Magnon interactions play a dual role. On the one 
hand, they renormalize magnon energy and, on the other 
hand, they can lead to magnon decay and result in a fi- 
nite lifetime. Both effects are treated on equal footing 
within the standard Green's function approach (Lifshitz 



and Pitaevskii 



1980 



Mahan 



2000}. The lowest-order 



self-energy diagrams produced by three- and four-particle 
vertices at T = are shown in Fig. [T] Analytic expres- 
sions for the diagrams in Figs. [TJa) and []Jb) are 



E a (k,cj) 
E 6 (k,a;) 



2 ^ UJ 

q 



n (1) (q;k)r 



■i0' 



(5) 



|^ 2, (q,k)| 



2 uj- 
q 



£k+q 



whereas the three-magnon decay process of Fig.JTJc) gives 
uj - e a - £ c 



E c (k,o;) = - > ! [ 



q-p 



(6) 



We included only independent momenta in the arguments 
of vertices in Eqs. ([5} and 

In the lowest order of perturbation theory, the real part 
of the on-shell self-energy E(k, uo = £\r) gives a correction 
to the dispersion 5e^ and the imaginary part yields the 
decay rate Tk = — Im£(k, £k). The source diagrams in 
Figs, [ljb) and [ljd) have no imaginary parts, while the 
magnon damping resulting from the decay diagram in 
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Fig. [TJa) is given by 

r k } = \ E l^3 (1) (q; k)| 2 % k - e q - e k - q ) • (7) 



The three-particle decay process in Fig.[TJc) yields a sim- 
ilar expression 



rL 4) = jEI^ (1) (^p; k )l 2 ^ 



Existence of the decay amplitudes and V± L) is a 
necessary, but by no means sufficient condition for spon- 
taneous magnon decays. The decay rate Tk is nonzero 
only if the energy conservation is satisfied for at least one 
decay channel. A general analysis of the kinematic con- 
ditions following from the energy conservation for two- 



^k-q-p) • 

(8) 



particle decays is given in Sec. Ill 



D. Three-boson interactions in spin liquids 

Apart from magnetically ordered antiferromagnets, 
there is a wide class of quantum-disordered magnets with 
spin-liquid-like ground states. At zero temperature such 
magnets remain completely isotropic under spin rotations 
and their ground-state wave function is a singlet state of 



the total spin £ t ot = 0; see Sec. |VII| for more detailed 
discussion and references. In spin-dimer systems and in 
some gapped chain materials, the low-energy excitations 
are 5 = 1 magnons separated by a finite gap from the 
singlet ground-state. This means that at any given mo- 
mentum k there is a triplet of excited states with the 
same energy e^. We use bosonic operators t\. a and t\^ a 
with a = x,y, and z to describe creation or annihila- 
tion of such quasiparticles, which are also called triplons. 
Having an intrinsic quantum number, spin polarization 
a, does not change appreciably the structure of the quar- 
tic terms in Eqs. and the SO(3)-invariant anhar- 
monic interactions are constructed by making all pos- 
sible convolutions for two pairs of polarization indices. 
An interesting observation is that in this case the spin- 
rotational symmetry does not forbid the cubic vertices. 
The total spin is conserved since two spin-1 magnons can 
form a state with the same total spin S tot = 1- The in- 
variant form of the cubic interaction for triplons is then 
given by 



%=|E ^(ki , k 2 ; k 3 ) e<^ (C*l a 0*ks7+ h - c -) 



(9) 



where the conservation of total spin is imposed by the 
antisymmetric tensor e a ^ 7 . In addition, there are certain 
lattice symmetries affecting the structure of the decay 
vertex in Eq. ([9|, which may still forbid the two-particle 
decays of triplons. Specific examples of that are discussed 
in Sec. ETD 



III. KINEMATICS OF TWO-MAG NON DECAYS 

The aim of this section is to consider kinematic con- 
straints in the two-particle decay process that follow from 
energy conservation: 

(10) 



For a magnon with the momentum k, the relation (10) 



is an equation for an unknown q with k being an ex- 
ternal parameter. Solutions of Eq. (10) form a decay 



surface, the locus of momenta of quasiparticles created 
in the decay. As a function of the initial momentum k 
the decay surface expands or shrinks and may disappear 
completely. In the latter case, the magnon with momen- 
tum k becomes stable. The region in k space with stable 
excitations is separated from the decay region where de- 
cays are allowed by the decay threshold boundary. 

For any given k the two-particle excitations form a 
continuum of states in a certain energy interval 

Ef a (k) < E 2 (k, q) = e q + e k _ q < £ 2 max (k) . (11) 

With the help of this statement it is straightforward to 
see that nontrivial solutions of Eq. ( 10 ) exist only if the 
single-particle branch and the two-particle continuum 
overlap. From this perspective it also becomes evident 
that the decay threshold boundary must be the surface 
of intersections of the single-particle branch with the 
bottom of the continuum Elf 1111 (Is). 

Suppose for a moment that spontaneous two-particle 
decays are completely forbidden, i.e., the one-magnon 
branch lies below the two-magnon continuum in the 
whole Brillouin zone 



for V k, q , 



(12) 



where the equality may be realized only for a trivial so- 
lution q = 0, provided that eo = 0. Applying the same 
relation (12) to £k- q on the right-hand side gives 

. for Vq,p, (13) 



£k < 



£k- 



which means that the one-magnon branch also lies below 
the three-magnon continuum as well as any n-magnon 
continua. Thus, if two-particle decays of Eq. (10) are 



prohibited for any value of k, the energy and momentum 
conservation also forbid decays in all n-particle channels 



with n > 3 (Harris et al, 1971). This is why it is im- 



portant to analyze kinematic conditions for two-particle 
decays as the first step even if the cubic terms are not 
present. 

Aside from finding whether spontaneous decays exist 
or not, there is another reason for considering kinematics 
of decays. In vestigati ng tw o-rot on decays in the super- 
fluid 4 He, Pitaevskii (1959) found that the enhanced den- 
sity of states (DoS) near the bottom of the two-particle 
continuum may produce strong singularities in the single- 
particle spectrum at the decay threshold boundary. Ad- 
ditional singularities may also occur inside the decay re- 
gion due to topological transitions o f the decay surface 
( |Chernyshev and Zhitomirsky| |2006[). A general analysis 
of these two effects is given in Sec. III.C 
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FIG. 2 (color online). Magnon decay threshold boundaries 
for the cubic antiferromagnet in an external field shown in 
one-eighth of the Brillouin zone for several values of the field. 
Decays are possible within the part of the Brillouin zone en- 
closing the point. 



A. Decay threshold boundary 



The resultant decay region is given by the union of all 
subregions enclosed by the boundaries obtained by con- 
sidering all the decay channels in (i)-(iv). Some of the 
calculated surfaces define the decay boundary while the 
rest correspond to saddle-points inside the decay region. 
Figure [2] shows an example of the decay threshold bound- 
aries for the cubic Heisenberg antiferromagnet in applied 
magnetic field, the problem discussed in more detail in 



Sec. [TV] Magnetic fields are given in units of the satura- 
tion field H s and below H* « 0.76H S decays are strictly 
forbidden. For H > H* the decay region grows out of the 
point, filling out the whole Brillouin zone at H = H s . 



B. Decay of acoustic excitations 

For a sufficiently complicated dispersion law the de- 
cay threshold boundary has to be found numerically by 
solving Eq. (10) and Eqs. (14)-(17). However, it is often 



possible to draw conclusions about the presence of decays 
by analyzing only the low-energy part of the spectrum. 
Such an approach is frequently use d in theoretical studies 
of phonon decays (iKosevich 20051) and we briefly review 



the corresponding arguments, which will also be relevant 
for the subsequent treatment of magnon decays. 



The bottom of the two-particle continuum can 
be found by imposing the extremum condition 
Vq£?2(k, q) = 0, which is satisfied if the velocities of 
two magnons are equal, 



v k - 



(14) 



Using this, the decay threshold boundary is determined 
by the system of equations (10) and (14). One can dis- 



tinguish several types of g eneral solutions, or decay chan- 
nels, for such boundaries (Chernyshev and Zhitomirsky 



2006| |Lifshitz and Pitaevskii[ |1980| ): 
(i) The emission of two equivalent magnons with equal 
momenta q = k — q = |(k + G), where G is a reciprocal 
lattice vector. In this case Eq. ( [l4| ) is satisfied automat- 
ically while Eq. ( 10 ) is reduced to 



(k+G)/2 



(15) 



(ii) The emission of an acoustic magnon with q — > 0. In 
this case Eq. (10) is fulfilled while Eq. (14) is equivalent 
to 



I vol 



(16) 



(iii) The emission of an acoustic magnon with q — >• 
(for magnets with several acoustic branches): 



£k = £k-Q; 



(17) 



(iv) The emission of two finite-energy magnons with dif- 
ferent energies but the same velocities. No simplification 
occurs in this case. 



1. Single acoustic branch 

In the decay process of a long-wavelength acoustic ex- 
citation the incident and the two emitted quasiparticles 
have the same velocity c. The linear approximation for 
the excitation energy, = ck, does not provide sufficient 
information on the possibility of decays and one has to 
include a nonlinear correction to the linear dispersion: 



£k ~ ck + Oik 



(18) 



For a weak nonlinearity ak 3 <C cfc, the momenta of de- 
caying and emitted quasiparticles are nearly parallel such 
that 



kqtp 2 



2(k-q) ' 



(19) 



where ip is a small angle between q and k. The energy 
conservation ( 10 ) within the same approximation is 



3akq(k — q) 



ckqcp 2 
2(fc - q) 



(20) 



The nontrivial solutions q, ip exist only for the posi- 
tive sign of the cubic nonlinearity a. While the asymp- 
totic form ( 18 ) is relevant to many physical examples, an 



even more general condition can be formulated to include 
other cases of gapless energy spectra. If the low-energy 
part of the spectrum is a concave function of the momen- 
tum d 2 e/dk 2 < 0, Fig. [3^ a), magnons remain stable. On 
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FIG. 3 (color online). Sketches of the quasiparticle spectra 
£k with (a) negative and (b) positive curvature of the acous- 
tic mode, respectively, (c) Spectrum with several acoustic 
branches. 



the other hand, for the upward curvature d 2 e/dk 2 > 0, 
Fig. |3|b), the low-energy excitations are unstable with 
respect to spontaneous decays. 

Note that according to the above criterion two-magnon 
decays are kinematically allowed in Heisenberg ferromag- 
nets. Indeed, the energy conservation equation (10) can 
be satisfied for quasiparticles with a ferromagnetic dis- 
persion law £k cx k 2 . However, as discussed in Sec. II, 



the decay vertices and V± L) vanish exactly in the 
Heisenberg ferromagnet and magnons remain well de- 
fined. Nevertheless, decays do exist in anisotropic planar 
ferromagnet s, for which the magnon numbe r conservation 
does not hold anymore and V} 1] is nonzero ( Stephanovich 
and Zhitomirskyl 120111 iVillainl 11974). 



(i) 



2. Several acoustic branches 

Lattice vibrations in crystals have three acoustic 
branches with different sound velocities corresponding to 
one longitudinal and two transverse phonons. Heisen- 
berg antiferromagnets with a helical or spiral magnetic 
structure also have three Goldstone modes, Fig. |3|c), at 
k = and ±Q, where Q is the ordering wave vector; see 
Sec. |V.A for more details. Energy conservation for the 
decay of a fast quasiparticle into two slow quasiparticles 
can be satisfied already in the linear approximation for 
the energy 



c\k = c 2 q + c 3 |k ■ 



for ci > c 2 ,cs 



(21) 



where all momenta are measured relative to the cor- 
responding Goldstone Q points. Thus, if this decay 
channel is compatible with momentum conservation, i.e., 
Qi = Q2 + Q3, then the fast excitation can decay re- 
gardless of the curvature of the corresponding acoustic 
branch. 



C. Decay singularities 

We now demonstrate the fact that the decay threshold 
boundary should generally correspond to a nonanalyt- 
icity in the particle's spectrum because of its coupling 
to the two-particle continuum. Consider the decay self- 
energy E a (k, uo) from Eq. ([5| in the vicinity of the thresh- 
old boundary for two-particle decays. Generally the de- 
cay vertex has no additional smallness in terms of the mo- 
menta of participating quasiparticles, i.e., vS^ = O(J). 



Focusing on the 2D case and expanding Eq. ( 10 ) in small 
Ak = k — k* and Aq = q — q* , where k* is the crossing 
point of the single-magnon branch and the bottom of the 
continuum and q* is the minimum point of i^(k*,q), 
energy conservation gives 

£k - £ q - £k- q « (vi - v 2 ) • Ak - ^ - ^ = , (22) 

where vi and w 2 are the velocities of the initial and final 
magnons and a and b are constants. Next, the singular 
part of the self-energy is given by 



£(k,e k ) cx J 



d 2 q 



fa - v 2 )Ak - q 2 x /a 2 - q%/b 2 + iO ' 



(23) 



and a straightforward integration in Eq. (|23| yields a 



logarithmi c singularity in the spectrum ( Zhitomirsky and 
Chernyshev[ |1999|) 



ReS(k,e k ) ~ In 



A 

~Ak\ 



r k = -ImS(k,e k )~9(Afc) : 



(24) 

where Q(x) is the step function. The imaginary part of 
E(k, £k) m Eq. ^ is directly related to the two-particle 
density of states. Therefore, in two dimensions it is natu- 
ral to have a jump in Tk upon entering the continuum, as 
obtained in Eq. ( 24 ) . In effect this demonstrates that the 



van Hove singularity in the continuum's density of states 
gets transferred directly onto the quasiparticle spectrum 
via the three-particle coupling. 

Consider now the other type of 2D van Hove singu- 
larity, a saddle-point inside the continuum. Expansion 
of the energy conservation in the vicinity of the saddle 
point yields the same form as in Eq. (22) with a change 
of sign in front of either Aq 2 or Aq 2 ^ Therefore, the 
singular part of the self-energy is 



£(k,e k ) cx 



d 2 q 



{v 1 -v 2 )Ak-ql/a 2 + ql/b 2 



iO 



(25) 



and the logarithmic singularity appears now in the imag- 
inary part: 



ReS(k,e k ) ~ sign(Afc) . 



r k ^ In 



A 

Afcl 



(26) 



again in agreement with the relation between the 2D two- 
magnon DoS and Tk- 
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Geomet ric consideration of the decay surface near sin- 
gularities (Chernyshev and Zhitomirsky , 2006) offers a 
useful alternative perspective: while crossing of the decay 
threshold boundary obviously corresponds to nucleation 
of the decay surface, the saddle-point singularity corre- 
sponds to splitting of the decay surface into two disjoint 
pieces. Altogether, both types of singularities (24) and 



( 26 ) result from topological transitions of the decay sur- 



face with a change in the number of connected sheets. 

In 3D the decay self-energy generally exhibits a less 
singular but still nonanalytic behavior. The logarith- 
mic peak at the decay threshold bounda ry is replaced 
by a continuous square-root behavior ~ y^\Ak\ for both 
ReS(k, £k) an d Tk- A remarkable exception is the thresh- 
old to the two-roton decay in the spectrum of superfluid 
4 He. The dispersion of roton excitations has a minimum 
on a sphere in the momentum space, leading to the same 
logarithmic divergence of the self -energy as in Eq. ([2 4]) at 
the energy twice the roton gap ( |Lifshitz and Pitaevskii| 
1980| ). This strong anomaly produces the "endpoint" in 
the dispersion cu rve £k of the superfluid 4 He predicted by 
Pitaevskii ( 1959 ). Such a termination point was later ob- 



served in inelastic neutron-scattering experiments ( Glyde 



et 



1998 Graf et al, 1974) 



The analysis of this section can be straightforwardly 
extended to the three-particle decay self-energy E c (k, oS) 
given by Eq. The general outcome is similar to the 
effects of higher dimensions: singularities are replaced 
by continuous albeit nonanalytic behavior. Overall, the 
two-magnon decays play a more significant role than the 
higher-order decay processes. Decays and associated sin- 
gularities are also enhanced by the low dimensionality of 
a magnetic system. 



IV. FIELD-INDUCED DECAYS 

We devote this section to spontaneous decays in the 
square-lattice Heisenberg antiferromagnet in strong ex- 
ternal field. In zero field the ground state of this model 
is a two-sublattice Neel structure, which, due to its 
collinearity, prohibits cubic anharmonicities and two- 
magnon decays; see Sec.[TTjB. The three-magnon and all 
other n-magnon decays are forbidden by energy conserva- 
tion, a s in the cubic-lattice model studied by |Harris et al. 



Ill 



( |1971| ) and in agreement with the arguments of Sec. 
Thus, in zero field and at T = magnons have an infi- 
nite lifetime. In the opposite limit, when spins are fully 
polarized by external field, the ground state is equivalent 
to that of a ferromagnet with all particle-nonconserving 
terms strictly forbidden so that the magnons are well de- 
fined again. It comes as a surprise that there exists a 
regime between these two limits where magnons become 
heavily damped throughout the Brillouin zone. 

Next we discuss general arguments for the existence 
of the field-induced decays in a broad class of antiferro- 
magnets and provide details of the spin-wave approach 
to the square-lattice antiferromagnet in a field. Focus- 




FIG. 4 (color online). Field- induced canted spin structure in 
the square-lattice Heisenberg antiferromagnet. 



ing mostly on the 5 = 1/2 case, we demonstrate that 
magnons become damped and even overdamped in most 
of the Bril louin zone for a range of fields near th e satura- 
tion field ( jZhitomirsky and Chernyshev , 1999). Confir- 
mation of these results from numerical (Syljuasen, 2008) 



and experimental studies is highlighted. We also illus- 
trate the occurrence of the threshold and saddle-point 



singularities discussed in Sec. |III| C in the perturbative 
treatment of the magnon spectrum and offer two ap- 
proaches to regularize them self-consistently. Detailed 
results of one such self-consistent regularization method 



(Mourigal et al. 2010) are presented and an example 



of the higher-dimensional extension of the square-lattice 
case is touched upon (Fuhrman et al., 2012). 



A. General discussion 

The evolution of an ordered magnetic structure of the 
two-sublattice antiferromagnet in an external field is rel- 
atively trivial. Spins orient themselves transverse to the 
field direction and tilt gradually towards its direction un- 
til they become fully aligned at H > H s , where H s is re- 
ferred to as the saturation field, see Fig.[4j In the absence 
of anisotropy and frustration there is no other phase tran- 
sition in the whole range from H = to H = H s . As first 



argued by Zhitomirsky and Chernyshev (1999), it is the 



dynamical properties of quantum antiferromagnets that 
undergo an unexpectedly dramatic transformation on the 
way to full saturation. The phenomenon responsible for 
this transformation is spontaneous magnon decay, which 
becomes possible in sufficiently strong magnetic fields. 

The field-induced magnon decays are generally present 
in a broad class of quantum antiferromagnets. Indeed, 
three-magnon interactions that couple one- and two- 
magnon states always exist in canted antiferromagnetic 
structures created by an external field; see Sec. [TT[ Usu- 
ally, kinematic constraints prevent spontaneous decays 
from taking place in zero or weak magnetic fields. It is 
easy to see that from the acoustic branch of the magnon 
spectrum (18), which retains the concave shape (a < 0) 



in small fields and thus is stable with respect to spon- 
taneous decays according to Sec. |III| B. However, in a 
strong enough magnetic field the convexity of the acous- 
tic spectrum must change. This is because at the satu- 
ration field, H = H s , the magnon velocity vanishes and 
the dispersion at low energies is that of a ferromagnet, 
£k oc /c 2 , which has an upward curvature. By continuity, 



the spectrum has to preserve its positive curvature for 
a certain range of magnetic fields H* < H < H s , where 
the threshold field H* corresponds to a = 0. This implies 
that the two-particle decays are present at H > H* . 



B. Square-lattice antiferromagnet 



The Heisenberg model on a square lattice, a toy model 



in the early days of the theory of aniferromagnets (Ander- 



son 



1952||Kubo||1952[ ), became a paradigmatic model in 



the late 1980s due to its direct relevance to the parental 



materials of the high-T c cuprates ( Chakravarty et al. 



1989 Manousakis 1991). It is also viewed as one of 
the prototypical strongly correlated models in condensed 
matter physics, used for benchmarking various theoret- 
ical meth ods against the more precise numer i cal tech- 
niques ( |Sandvik[ [l997| |White and Chernyshev] |2007 ). 



The standard spin-wave theory provides an accurate 
description of the static and dynamical properties of 
the square-lattice antiferromagnet in zero field even for 



5=1/2 flHamer 


et al] 1992 Igarashi 




1992; Igarashi 


and Nagao[ |2005 




Weihong and Hamer 




1993|. In the 


words oflSandvik 




1997): "The presently most accurate 



calculations [for the 5 = 1/2 case]... indicate that the 
true values of the ground-state parameters deviate from 
their 1/S 2 spin- wave values by only 1-2% or less." The 
spin- wave calculati ons also offer a fine fit to the overall 
magnon dispersion ( Igarashi and Nagao, 2005; Syromyat 



nikov 2010 Weihong and Hamer[ 1993 ) with the only de- 



viatio n from the numerical results along the (7r,0)- (0, tt) 



path (Sandvik and Singh 2001 Zheng et al. 



2005) 



Heisenberg models in strong fields have received some- 



what less attention u ntil recently [see, however, |de Jongh 
|and Miedema| ( |1974D and |de Groot and de Jongn|Q1986| )], 



primarily because the majority of available materials had 
exchange constants much larger than the strength of 
available magnetic fields. The constraints are even more 
stringent on neutron-scattering studies in the field, where 
the practical limit is currently at about 14T. Neverthe- 
less, recent developments in the synthesis of molecular 
based antiferromagnets with moderate strength of ex- 



change coupling between spins (Coomer et al. , 2007 Lan- 
caster et al] fZU07\ |Woodward et a/.| |2002| |Xiao et al] 



2009 ) have opened the high-magnetic-field regime to ex- 
perimental investi gations for a number of layer ed square- 
lattice materials (iTsyrulin et al] 120091 120101). Further- 



more, new field-induced dynamical effects can be present 
in antiferrom agnets with other lattice geometries ( |Coldea 



et al, 2001). 



Another class of antiferromagnets directly relevant to 
our discussion incorporates quantum spin systems with 
singlet ground states and gapped triplet excitations. 
They are often called BEC magnets because the field- 
induced ordering in them can be described in terms of 
the Bose-Einstein condensation (BEC) of triplet excita- 
tions. In recent years these quantum spin-gap magnets 
have been intensively studied both experimentally ( |Jaime| 



et al, 2004 Nikuni et al 


2000 


Regnault et al, 2006) 


and theoretically flAmeck 


1991 


Affleck and Wellman| 


1992 Giamarchi and Tsvel 


ik, 1999; Mila, 1998|) and con- 


tinue to attract attention (|Giamarchi et al] [2008 ) . We 



decays discussed here equally applies to the vicinity of 
the triplet condensation field H c in these magnets be- 
cause of the duality between H c and the saturation field 
H s (Mila 1998). This is highly advantageous because in 
many BEC magnets the lower condensation field is read- 
ily accessible, which makes them prime candidates for 
observing magnon decays directly in inelastic neutron- 
scattering experiments. It can be argued that the re- 
cently observed suppression of thermal conductivity in 



the vicinity of critical fields in one such material ( Ko- 



hama et al. 2011| ) is due to magnon decays. 



Extension of the spin-wave theory to finite magnetic 
fields was develop e d by Osano et al.\( 1982 ^_ Zhitomirsky| 
and Nikuni ( 1998 ), Zhitomirsky and C herny shev| fll 999 ) , 



Chernyshev and Zhitomirsky (2009b), and Mourigal et 



al. (2010). For alternative approaches applicable in the 



vicinity of the saturation field, see also Batyev and Bra 



ginskii (1984), |Gluzman| (|1993[), Kreisel et al. (2008), 



and Syromyatnikov| ( |2009| ). Close agreement of the spin- 



wave calculations with exact diagonalization and quan- 
tum Monte Carlo (QMC) results for the ground-state 
properties in the w hole range of fields < H < H s 
was established by Liischer and Lauchli (2009). The 



same work has also confirmed the earlier QMC results 
by Syljuasen ( 2008[ ) that demonstrated clear signatures 
of magnon decays in the high- field regime. We provide 
more details on that in the following. 



1. Model and formalism 

We begin with the Heisenberg Hamiltonian of nearest- 
neighbor spins on a square lattice in a magnetic field 
directed along zq axis in the laboratory reference frame, 



U 



(ij) 



(27) 



Here, J is the nearest-neighbor coupling and H is the 
field in units of g\iB- With the details of the technical 
appr oach exp licated by Zhitomirsky and Nikuni (1998) 
and Mourigal et al. ( 2010| ), we summarize here the key 
steps of the spin-wave theory approach to this problem. 

First we align the local spin quantization axis on each 
site in the direction given by the canted spin configu- 
ration shown in Fig. [4j The corresponding transforma- 
tion of the spin components from the laboratory frame 
(#o, 2/o> z o) to the local frame is 

Sf = Sf sin + Sf e iQ ' r * cos , 5f = Sf , 

Sf = -Sf e iQ r * cos + Sf sin , (28) 

where is the canting angle to be defined from the energy 
minimization and Q = (7r, tt) is the ordering wave vector 
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of the two-sublattice Neel structure. Next the standard 
Holstein-PrimakofT transformation bosonizes the spin op- 
erators Sf = S — a\a,i, S~ — a\(2S — aja^) 1 / 2 , with 
Sf = Sf±iSf. Expanding square-roots, one obtains the 
bosonic Hamiltonian as a sum, 



(29) 



each term being of the nth order in bosonic operators 
a] and a$ and carrying an explicit factor S 2 ' 71 ^. This 
form provides the basis for the 1/5 expansion. Important 
consequences of the noncollinear spin structure are the 
terms SjS^ y \ which describe coupling between trans- 
verse S x ( y ^ and longitudinal S z oscillations and yield cu- 
bic anharmonicities (H3) in the bosonic Hamiltonian. 



Minimization of the classical energy 7-Lq in Eq. (29) 
fixes the canting angle of the classical spin structure to 
sin (9 = H/H S1 where the saturation field is H s = 8JS. 
This procedure also cancels the Hi term that is linear in 
a\ and a^. After subsequent Fourier transformation, the 
harmonic part of the Hamiltonian (H2) is diagonalized 



by the Bogolyubov transformation, = u^b^ 
yielding 



v k bl 



k' 



Usw = Y, i * h lK ( 3 °) 

k 

+ 1 E v s {1) ( k > <0 ( b i-^Q h iK + hx.) + . . . , 

k,q 

where = + fek with being the magnon dis- 
persion given by the linear spin-wave theory and fek 
is from the 1/5 corrections due to angle renormaliza- 
tion and Hartree-Fock decoupling of cubic and quartic 
terms in Eq. p9| ), respectively. The ellipsis stands for 

(2) 

the classical energy, the three-boson source term V3 as 
in Eq. ([3|, and the higher-order terms in the 1/5 ex- 
pansion. Although some of these terms do contribute to 
the subsequent results, this abbreviated form constitutes 
the essential part of the Hamiltonian and is sufficient for 
our discussion. For the e xplicit expressions of S e^ and 



vertices Vo 1 ^ and V** } , see 



(2) 



2. Decay regions 



Mourigal et al (2010) 



One can already gain detailed insight into the decay 
conditions in the manner discussed in Sec. |III| by analyz- 
ing the magnon spectrum in the harmonic approxima- 
tion, which is given by 



£k = 4J5V(l + 7k)(l-cos2# 7k ) 



(31) 



where 7k = |(cos /c^ + cos k y ). The acoustic mode near Q 
follows the asymptotic form (18) £Q+k ~ ck + ak 3 with 



2V2JS cos (9, 



c ( 9 „ 9 + cos4(^\ . 
16 H ° 6-^>' (32) 
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FIG. 5 (color online). Decay regions in one- fourth of the 
Brillouin zone for representative values of H/H s indicated by 
numbers. Solid and dashed lines are threshold boundaries for 
two- and three-magnon decays, respectively. 



where ip is the azimuthal angle of k = (k Xl k y ). The 
curvature of the spectrum a changes its sign for k along 
the diagonal (cp = 7r/4) at the decay threshold field 



H* 



V7 



H s « 0.7559 H 8 



(33) 



Note that this expression remains valid for the cubic- 
and for the layered square-lattice antiferromagnet with 
arbitrary antiferromagnetic interlayer coupling. 



Staggered canting of spins in Eq. (28) "shifts" the mo- 
mentum in the two- magnon decay condition in Eq. (10), 
which now reads as = £ q + £k-q+Q- Using the Bose- 
condensate perspective of Sec. |II.B| such a change is nat- 
ural as the magnon going into the condensate now carries 
the momentum Q; see Eq. Q. Apart from that the kine- 
matic analysis of Sec. |III| remains the same. Numerical 
results for the decay region at H > H* are shown in 
Fig. [5] Up to H « 0.84i^ s the boundary of the decay 
region is determined entirely by the decays into a pair 
of magnons with equal momenta; see Sec. |III.A[ type 
(i) in Eq. ([15]). At fields higher than 0.85i^ s the decay 



channel with emission of an acoustic magnon, type (ii) 
in Eq. (16), creates protrusions leading to more compli- 



cated shapes of the decay region. Above 0.9i^ s the decay 
region spreads almost over the entire Brillouin zone, sim- 
ilarly to the 3D decay regions in Fig. [2] In Fig. [5j we also 
show boundaries for the three-magnon decays, which are 
less restrictive and cover larger areas of the Brillouin zone 
than their two-magnon counterparts. 



C. Spectral function and self-consistent Born 
approximation 

Information on the magnon energy renormalization 
and decay-induced lifetime due to interactions can 
be obtained from the single-particle spectral function 
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FIG. 6 (color online). Magnon spectral function for the spin- 
1 antiferromagnet at H = 0.9H S for k = 7r(g, q). Dashed and 
solid lines correspond to Born approximation and restricted 
self-consistent Born approximation, respectively. From [Zhit-| 
omirsky and Chernyshev ( 1999| ). 



A(k, lj) = — (1/tt) ImG(k, lj). Here G(k, lj) is the inter- 
acting Green's function, 



G 1 (k,o;) = cj -e k - E a (k,cj) - E 6 (k,cj) 



(34) 



where E a? b(k, a;) are the decay and source self-energies 
shown in Fig.jl] with the "one-loop" expressions for them 
given in Eq. (pi. Importantly, the magnon spectral func- 
tion is directly related to the dynamical structure factor 
ff(k, lj) measured in n eutron-scattering experiments; see 
Mourigal et al (2010) for the details of their connection. 



The self-energies in Eq. (34) originate from coupling 



to the two-magnon continuum and already in the lowest 
order one expects the spectral function to exhibit an inco- 
herent component due to two-magnon states in addition 
to the quasiparticle peak that is broadened due to de- 
cays. The Born approximation (BA) results for A(k, lj) 
for several momenta along the TM line at H = 0.9 H s 
are shown in Fig. [6] by dashed lines. While A(k,u) ex- 
hibits incoherent subbands, the peaks survive even for 
k deep inside the decay region in Fig. [U This is due 
to the lack of self-consistency in the lowest-order BA, as 
one- and two-magnon states are treated with different 
accuracy. Specifically, the magnon energy is renormal- 
ized downward by the real part of E a? 5(k,cj) in Eq. (34), 



whereas the energy of the continuum is still calculated 
in the harmonic approximation. Taking into account en- 
ergy renormalization of the continuum means "dressing" 
of the inner lines of the decay diagram in Fig. [I] The 
main technical difficulty with such a self-consistent Born 
approximation (SCBA) is that it opens an unphysical gap 
for long-wavelength magnons in violation of the Gold- 
stone theorem. 

This difficulty notwithstanding, the restricted SCBA 
scheme was implemented by |Zhitomirsky and Cherny-| 
shev (1999), which included partial renormalization of 
the two-magnon continuum. It enforced the acoustic 



form of the long- wavelength magnon dispersion (Jl8| and 
incorporated the self-consistent renormalization of one 
inner magnon line in the decay diagram of Fig.[TJa). The 
results of the numerical solution of Dyson's equation in 
the restricted SCBA are shown in Fig. [6] by solid lines. 
Not only are the peaks that were artifacts of the non- 
self-consistent treatment completely wiped out, but the 
overall values of magnon damping are large away from 
the Q and T points, the latter corresponding to the uni- 
form precession mode which is stable. 

These results received strong confirmation from the nu- 



merical QMC simulations by Syljuasen (2008) (see Fig. [71 
and, subs equently, from the exa c t diag onalization (ED) 
study by |Liischer and Lauchli| ( |2009[ ). Both numeri- 
cal works found significant modifications in the magnon 
spectra above the threshold field: an explicit broaden- 
ing of peaks in the dynamical structure factor revealed 
by QMC data and a qualitative transformation in the 
structure of energy levels and in distribution of spectral 
weights observed in ED. Figure [7] shows the intensity plot 
of the S zz (q_, uj) component of the dynamical structure 
factor from Syljuasen] ( |2008[ ). Since S zz (q_, uj) is propor- 
tional to the spectral function with a shifted momentum 
A(q— Q,cj) (Fuhrman et al 2012), gapless and uniform 
precession modes trade their places in Fig. [7] 

There are close qualitative and quantitative similarities 
between the analytical and numerical results of Figs. [6] 
and [7[ Magnon decays and spectrum broadening are 
strongly pronounced and follow decay regions suggested 
by the kinematic analysis of Sec. IV. B. Syljuasen also em- 
phasized the nontrivial redistribution of spectral weight 
resulting in non-Lorentzian "double-peak" features in the 
dynamical structure factor, also seen in the analytical 
results in Fig. [6j albeit less pronounced. In spite of the 
approximations of the analytical approach and of the un- 
certainties in the QMC results due to numerical interpo- 




(1,0) (1,1) (1/2,1/2) 

q (jt/a,jt/a) 

FIG. 7 (color online). QMC results for the S zz (q, lj) compo- 
nent of the dynamical structure factor for the S — 1/2 square- 
lattice antiferromagnet at H — 0.875if s along the YXMY 
path. Black li nes are th e magnon linewidths in the Born ap- 
proximation \ Sylju asen 2008"). The box marks the k-cj region 
of Fig. [6] for comparison. From Olav Syljuasen, 2008. 
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FIG. 8 (color online). Decay region and associated singular- 
ities for H = 0.9H S . Solid, dashed, and dotted lines denote 
the decay channels into a pair of magnons with equal and 
different momenta, and emission of an acoustic magnon, re- 
spectively. Thick portions of solid and dashed lines indicate 
the saddle-point singularity in the continuum, while the thin 
parts correspond to local minima. 



lation from imaginary to real frequencies, such a corre- 
spondence between two very different methods is truly 
remarkable. 

Following the original prediction of the field-induced 
spontaneous magnon decays, there is an ongoi ng search 



for suitable spin-^ square-lattice materials (|Coomer et 



q/4 12007| [Tancaster et al\ |2QQ7| [Tsyrulin et a/.[|2UU9|) to 

inve stigat e corresponding effects! Recently, |Masuda e£| 
aT ( |2Q1Q[ ) found indications of the field-induced decays 



in the inelastic neutron-scattering spectra of the spin- 
layered square-lattice antiferromagnet Ba2MnGe207 



D. Decay singularities and SCBA 

While the preceding consideration demonstrates the 
spectacular role of spontaneous decays in reshaping the 
magnon spectrum in the high-field regime of the S = 
1/2 square-lattice antiferromagnet, several aspects of the 
bigger picture remain unclear. In particular, what about 
the threshold and saddle-point singularities discussed in 



Sec. III? Are there termination points in the magnon 
spectrum? How strong is the magnon damping in more 
quasiclassical, S > 1 , systems? Th e se qu estions were 
addressed recently by Mourigal et al. (2010). 



First we substantiate the generic analysis of the thresh- 
old singularities provided in Sec. |III| by giving a spe- 
cific example from the square-lattice antiferromagnet. In 
Fig. [8] we show the decay region for H = 0.9 H s to- 
gether with the set of contours obtained by the analysis 
of Sec. |III.A| The solid, dashed, and dotted lines are the 
singularity contours for the decay into magnons with the 
same (15) and different momenta, and for the emission 



(saddle point) of the continuum. The same singularity 
contour may correspond to a minimum in one part of the 
Brillouin zone and to a saddle point in the other. This 
determines whether the logarithmic singularity will oc- 
cur in the real or imaginary part of the self-energy, as 
discussed in Sec. IIII.CI 

Within the spin- wave expansion, the 1/S renormaliza- 
tion for the magnon spectrum is given by 



£k = £k + Re E a (k, e k ) + E 6 (k, e k ) 



(35) 



where £k = £k + fek already contains the Hartree- 
Fock corrections to the linear spin-wave dispersion £ k 
[see Eq. (30)], and the magnon damping is T k = 
-ImS a (k,£k) in the same approximation. In Fig. [9] we 
present numerical results for the renormalized magnon 
energy and damping for H = 0.9i7 s , S = 1/2, and for 
k along the MTX line, which crosses several threshold 
boundaries according to Fig. [8j The magnon energy ex- 
hibits singularities that are clearly related to spontaneous 
magnon decays: every spike or jump in £ k is accompanied 
by a jump or a peak in the decay rate T k . In agreement 
with Sec. |III.C[ the logarithmic singularities occur on a 
crossing of the singularity contour with substantial den- 
sity of two-magnon states, a minimum or a saddle point 
of the continuum outlined in Fig. [§J 

Such singularities invalidate the usual 1/S expansion 
because the renormalized spectrum contains divergences. 
The encountered problem is generic and must be common 
to a variety of the 2D models regardless of the value of the 
on-site spin S. For systems with large spin the situation 
is especially aggravating as the 1/S expansion changes 
from a reliable approach to one producing unphysical di- 
vergences, undermining any sensible comparison with ex- 
perimental or numerical data. 

Singular corrections to the spectrum obtained in the 
leading order require a higher-order regularization, which 
depends on the form of e k . For example, two-roton decay 
in the superfluid 4 He yields logarithmic divergence of the 
one-loop diagram in Fig. [TJa) near the decay threshold. 
Using regularization of the two-particle scattering ampli- 




of an acoustic magnon (16), respectively. Line thickness 



indicates a singularity that corresponds to a minimum 



FIG. 9 (color online). Magnon energy and decay rate in units 
of J for the spin-| square- lattice antiferromagnet along the 
MTX path at H = 0.9H S . The two upper curves show the 
harmonic £ k , dashed line, and the renormalized magnon en- 
ergy £k, solid line. The lower curve gives the decay rate r k . 
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FIG. 10 (color online). Intensity maps of the magnon decay 
rate Tk calculated within the iSCBA for several values of ex- 
ternal magnetic field and for two values of the spin: S = 1 
(upper panels) and S = 5/2 (lower panels). 



tude, Pitaevskii (1959) showed that such thresholds be- 



come termination points of the spectrum; this was later 
confirmed experimentally ( |Glyde et al. , 1998). However, 
this salient feature of the 4 He spectrum is owing to the 
form of £k? which ensures that rotons created in the decay 
process are well defined. The problem considered here is 
qualitatively different because magnons created in an el- 
ementary decay are also unstable. This can be seen from 
the evolution of the decay regions versus field in Fig. |5j 
which begins with the acoustic branch where the typi- 
cal decay processes produce magnons further within that 
region. 

Such "cascade" decays imply the following. 
First, the threshold and saddle-point singularities 
in Eqs. ( [24] ) and ([26]) are removed by changing 



\Ak\ 



v/(A&) 2 + (f q */c q *) 2 , where T 



q* 



the 



and 



damping of the "singular" d ecay products emu c q * 
is related to their velocities ( |Mourigal et al. , 2010| . 
Second, the higher-order regularization treatment must 
include damping of the decay products. Thus, as in 
the Sec. IV. C, dressing of the bare propagators in 
Fig. [TJa) within the SCBA should yield a regularized 
spectrum. Having in mind the technical difficulty of this 
approach we used instead another scheme, referred to as 
imaginary-part SCBA (iSCBA). This treats the magnon 
decay rate Tk self-consistentl y, but neglects q u antum 
corrections to £k altogether (Mourigal et al. 



2010). 



Indeed, away from singularities and for larger spins, 
S > 1, renormalization of the spectrum Se-^ ~ O(J) is 
small compared to the bare energy ~ O(SJ) and can 
be neglected. Although the intrinsic damping is of the 
same order Tk ~ O(J), it represents the most important 
qualitatively new feature of the spectrum acquired due 
to spontaneous decays. 

In the simplest realization of the iSCBA we impose the 
following ansatz on the Green's function: 



This approximation is equivalent to assuming a 
Lorentzian shape of the quasiparticle peaks in the dy- 
namical response. Neglecting contributions from the 
source diagram, the self-consistent equation for the 
magnon decay rate Tk becomes 



r k 



q 



|y 3 w (q;k)nr q + r k _ q+Q ) 
Ok - £ q - £k-q+Q) 2 + (r q + r k _ q+Q ) 2 



(37) 

Such a scheme was previ ously applied to the proble m of 
phonon damping in 4 He (Sluckin and Bowley 1974) and 
in quasicrystals ( |Kats and Muratov| 2005). 

One of the qualitative results of the iSCBA approach 
can be immediately anticipated: since by dressing the 
diagram in Fig. [TJa) we implicitly take into account all 
n-particle decay processes, the k-region where magnons 
are damped is no longer confined to the two-magnon de- 
cay region. In Fig. [5] we provide an indication of that 
by showing the decay regions for the three-magnon pro- 
cesses (dashed lines). These are generally weaker than 
the leading two-magnon decays, but their kinematics is 
less restrictive as for any given H > H* the two-magnon 
decay area in Fig. [5] is enclosed by the three-magnon one. 
In the limit of n — >> oo, the decay region extends farther 
and can cover the entire Brillouin zone in the high-field 
regime. 

The same conclusion is reached if we consider broad- 
ening of the decay threshold singularity by the finite life- 
time of the final- state magnons. This broadening implies 
that the boundary in the momentum space between sta- 
ble and damped magnons is smeared out and the decay 
probability disappears gradually rather than in a steplike 
fashion. 



The numerical solution of the iSCBA Eq. (37) for two 
representative values of spin 5 = 1 and S = 5/2 and for 



several field values is presented in Fig. 10 in the form of 
intensity maps of rk/4J5 vs k in the same part of the 
Brillouin zone as in Fig. [8] As discussed, the singularities 
in Tk are washed out and finite damping develops outside 
of the original decay regions of Figs. [5] and [8] Comparison 
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(36) 



of the upper (5 = 1) and lower (S = 5/2) rows in Fig. 
reveals similar patterns in the momentum distribution of 

r k . 

Specifically, in lower fields, H < 0.9i^ s , magnon damp- 
ing is most significant in the broad region around the 
(7r/2,7r/2) point, which is a consequence of the large 
phase-space volume for the two-particle decays in this 
region. In higher fields, the maxima in the decay rate 
shift to the TX (TX f ) lines where the enhancement of 
damping is due to the remnants of the strong logarith- 
mic saddle-point singularities in Fig. [9] For S ^> 1 one 
can show that Tk ~ 0{J) away from the former singu- 
lar contours but there is a parametric enhancement of 
damping Tk ~ OUjnS ) in their vicinity ( |Chernyshev 
and Zhitomirsky|[2009a[ ). 

To summarize, the self-consistently calculated magnon 
damping in the large-5 square-lattice antiferromagnet in 
strong field is free from singularities and does not ex- 
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FIG. 11 (color online). Intensity plots of the dynamical structure factor S(k,uA, for fields 0.9H S and 0.95ij" s , S = 1/2, and for 
the k path in the k z — plane, along the T— )>M— )>X— ^X ; — )>r line; see Fig.[2]for notations. From Fuhrman et al. (2012). 



ceed rj^ ax ~ 0.7 J — 0.8 J. In other words, spin waves in 
large-^ systems acquire finite lifetimes at H > H* but 
remain well-defined quasiparticles. The regularized sin- 
gularities lead to a parametric enhancement of magnon 
damping and produce distinct field dependence and mo- 
mentum dependence of T^. The above results are in qual- 
itative agreement with the inelastic neutron-scattering 
data on the spin-| layered square-lattice antiferromag- 



net Ba2MnGe207 (Masuda et al., 2010), although more 
refined measurements are necessary for a detailed com- 
parison between theory and experiment. 



E. Decays in three dimensions 

In the two preceding sections we presented results of 
self-consistent schemes that dealt with the unrealistic es- 
cape of the quasiparticle peaks from the unrenormal- 
ized two-magnon continuum and with singularities en- 
countered in the 1/5 expansion. However, as noted in 
Sec. III.C[ the decay-induced divergencies in Tk and £k 
in 2D are replaced by a contin uous albeit nonanalytic 
behavior in 3D. Recent work by Fuhrman et al. (2012) 
demonstrated that a moderate interlayer coupling mit- 
igates 2D singularities and represents an alternative to 
the above self-consistency schemes. This approach al- 
lows one to obtain the excitation spectrum using the 
standard spin-wave expansion, without the use of self- 
consistency, with no restrictions on the spectral shapes 
and at a fraction of the computational cost. As a specific 
example, we calculated the dynamical structure factor 
of the quasi-2D S =1/2 square-lattice antiferromagnet 
with moderate interlayer coupling J 1 j J — 0.2, which is 
relevant to one of the prospective materials for observing 
magnon decays, (5C AP)2CuCU, where 5C AP is 2-amino- 
5-chloropyridinium ( jCoomer et al. , 2007). 

In Fig. [TTJ we show these results in the form of inten- 
sity plots of the dynamical structure factor, 5(k, oj) for 
H = 0.9H S and 0.95i^ s along a representative path in 
the Brillouin zone; see Fig. [2] for notations and Fuhrman 



et al. ( |2012 ) for details. In addition to the main fea- 
tures that are directly related to the magnon spectral 
function A(k,u;), the out-of-plane transverse mode con- 
tribution is visible as a shadow, shifted by the ordering 
vector Q = (7r,7r,7r). We performed a Gaussian convolu- 
tion in the uj direction with a a width of 0.1 J. This step 
is intended to mimic the effect of a realistic experimen- 
tal resolution, from which we can also draw quantitative 
predictions of the relative strength of the quasiparticle 
and the incoherent part of the spectrum. 

Complex spectral lineshapes, very much distinct from 
the conventional quasiparticle peaks, are clearly visible 
in Fig. [TT] This demonstrates that the effects of spec- 
tral weight redistribution and broadening due to sponta- 
neous decays are substantial and should be readily ob- 
served in experiment. The structure factor exhibits a va- 
riety of unusual features, including the double-peak line- 
shape for the TM direction discussed previously. Along 
some of the k directions for H = 0.9H S , the renormal- 
ized quasiparticle peaks escape the unrenormalized two- 
magnon continuum and survive despite being formally 
within the decay region shown in Fig. [2] Yet these peaks 
are accompanied by continuumlike subbands that accu- 
mulate significant weight. Upon increase of the field, the 
continuum overtakes the single-particle branches in large 
regions of the Brillouin zone, washing them out and cre- 
ating more double-peak structures, e.g., very distinctly 
along the XT direction. There is a particularly spec- 
tacular advance of the continuum along the XX' line, 
where magnon broadening is most intense with dramatic 
spectral weight transfer to higher energies. 

The overall landscape of the high-field dynamical 
structure factor shown here is diverse and intriguing. 
Our results, consistent with prior numerical and ana- 
lytical studies, provide a comprehensive illustration of 
its details. Inelastic neutron-scattering measurements on 
suitable quantum antiferromagnets will be important to 
elucidate the accuracy of our theoretical results. 
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V. ZERO-FIELD DECAYS 

In this section we discuss spontaneous decays in non- 
collinear antiferromagnets in zero magnetic field. As ex- 
plained in Sec. [TTJ the noncollinear magnetic structures 
necessarily lead to the coupling of transverse and longi- 
tudinal spin fluctuations that translates into cubic an- 
harmonicities in the effective magnon Hamiltonian. This 
opens the door to spontaneous two-magnon decays with- 
out the assistance of an external field, provided the kine- 
matic conditions of Sec. Hill are fulfilled. Here we concen- 
trate on the case of spiral Heisenberg antiferromagnets, 
in which the noncollinear ground state spontaneously 
breaks the SO (3) symmetry of the spin Hamiltonian. 
We show that the kinematic conditions that allow two- 
magnon decays are always fulfilled in such a case, ow- 
ing to the existence of three acoustic branches in the 
excitation spectra. Next we elucidate these features and 
consider magnon decays in the spin-^ triangular-lattice 
antiferromagnet in more detail. 

Our subsequent discussion does not include non- 
collinear magnetic structures induced by the anisotropic 
Dzyaloshinskii-Moriya interactions, which have attracted 
a lot of interest because of experimental observation of 
the skyrmion lattices (RoBler et al. 2006). We also 



leave aside the other cases of noncollinear antiferromag- 
net s, such as kagome-lattice potassium jarosite ( |Matan 
al. , 2006[ ) and anisotropic triangular-lattice antiferromag- 
nets, discussed by |Chernyshev and Zhitomirsky ( |2009a ), 
although many qualitative aspects of our consideration 
below remain valid for them. Yet another class of zero- 
field spontaneous decays that do not require noncollinear- 
ity, the three-magnon decays in anisotropic ferromagnets 
dStephanovich and Zhitomirsky) |2011| |Villain| |1974[ ), is 
also left out. 



A. General discussion 

Noncollinear spin structures in isotropic antiferromag- 
nets appear as a consequence of magnetic frustration: 
either because of the geometry of the nearest-neighbor 
bonds as, e.g., for a triangular lattice, or due to compet- 
ing exchange interactions beyond the nearest-neighbor 
shell. Here we consider a general case of the Heisenberg 
Hamiltonian 



(ij) 



(38) 



where spins occupy the sites of a Bravais lattice and r = 
Yj — Yi. The lowest-energy classical state of this model 
is a spin helix with a wave vector Q, which is found by 
minimizing the Fourier transform of the exchange energy 
J q = E r J ( r ) cos(q-r). 

The spin-wave theory has been applied to spiral mag- 
netic structures in a number of works; see, for exam- 
ple, [Chubuk^ ([1984]); [Ohyama^ Simi- 
larly to the procedure discussed in Sec. IV. B the classical 



spin configuration defines the rotating local basis for spin 
quantization, which is related to the laboratory frame via 



Sf° = Sf sinQ-r* + Sf cosQ-r* , Sf° 
Sf = Sf cosQ • Yi - Sf sinQ • . 



SV 



(39) 



The subsequent 1/5 expansion follows standard steps 
(see Sec. IV. B[ ) and yields the spin-wave dispersion in 
the harmonic approximation: 



£ k = SyVk - J Q ) (hJQ+k + § Jq-k - Jq) • (40) 

One can see that for an arbitrary form of J q the exci- 
tation energy vanishes at three wave vectors: k = and 
±Q. The existence of these three Goldstone modes fol- 
lows from the complete breaking of the rot ational SO (3) 
symmetry in the noncollinear ground state ( Andreev and 



Marchenko 



1980; Dombre and Read , [1989 ) . For a planar 



magnetic structure, spin waves with k = ±Q correspond 
to the out-of-plane oscillation modes related by symme- 
try, whereas the k = mode represents the in-plane oscil- 
lations of the spin spiral. Overall there are three acoustic 
branches, shown schematically in Fig. [3jc), with two dis- 
tinct spin- wave velocities Co and cq. 

As was noted, coupling between one- and two-magnon 
states is inevitable in noncollinear structures. Hence, 
to verify the presence or absence of spontaneous decays 
in spiral antiferromagnets one only needs to analyze the 
kinematic constraints following from the energy conser- 
vation (10). Our discussion in Sec. III.B suggests that 



in a system with multiple acoustic modes the fast low- 
energy magnons can always decay into slow ones. Using 
Eq. (40) one can obtain the velocities of the two acoustic 



modes as 



c = SJA(J - J Q ) , A=- V^J, 



klk=Q 



cq = Sy A (| Jq + \ J2Q — Jq) . 



(41) 



Since the ordering wave vector of a spiral structure has 
to satisfy the inequality 2Q 7^ G^, the G^'s being recip- 
rocal lattice vectors, the two velocities in Eq. (41) are 



generally different. If Jo > J2Q, the in-plane mode is 
faster: Co > cq. For example, this relation holds if all 
exchange constants are antiferromagnetic, J(r) > 0. In 
such a case the decay channel k ^ (Q + q) + (— Q + q') 
is always open for sufficiently small k. Since for the lin- 
earized spectra decays would be permitted for any k, it 
is the short-range behavior of J q that limits the size of 
the decay region, which typically occupies an extended 
region of the Brillouin zone including the T point. Gener- 
ally, acoustic branches emanating from the ±Q points re- 
main undamped for an incommensurate Q, although, see 



the discussion of the triangular-lattice case in Sec. V.B 



If, on the other hand, Jo < J2Q, which happens when 
some bonds are ferromagnetic, out-of-plane modes are 
faster: cq < cq. Now the decay channel is switched to 
(Q + k) — » (Q + q) + q' for small k such that the decay 
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FIG. 12 (color online). The ordered 120° spin structure of 
the triangular-lattice Heisenberg antiferromagnet. 



region surrounds the ±Q points, while magnons in the 
center of the Brillouin zone remain undamped. 

In this discussion we disregarded possible anisotropy 
of the spin- wave velocities. However, since both cq and 
cq have the same angular dependence via A in Eq. (41), 
the anisotropy plays only a secondary role, which does 
not change the qualitative conclusion about the magnon 
decay region. 

A finite lifetime of magnetic excitations was noted in 
several studies of noncollinear antiferromagnets within 



the higher-order spin- wave theory ( Dalidovich et al 
20061 IQhyama and Shibaj [1993| |Starykh et al.\ |2QQ6 



Veillette et al. 



2005). The generality of the phe- 
nomenon of spontaneous magnon decays in spiral anti- 
ferromagnets was emp hasized in our work (Chernyshev 
|and Zhitomirs ky, 2006), in which we studied in detail the 
triangular-lattice Heisenberg antiferromagnet considered 
next. 



B. Triangular-lattice antiferromagnet 

The semiclassical S ^> 1 triangular-lattice antiferro- 



magnet orders in the 120° spin structure shown in Fig. 12 
For the case of S = 1/2, it was proposed that the en- 
hanced quantum fluctuations may completely destroy 



the long-range antiferromagnetic order (Anderson, 1973 



Fazekas and Anderson[ |1974[ ). However, the spin- wave 
calculations suggested that the 120° state is stable for 
S = 1/ 2, albeit with a significantly reduced ordered mo - 



ment (Jolicoeur and Le Guillou, 1989 Miyake, 1992). 



Th e exac t diagonalization of small clusters ( Bernu et 
al. 1994) together with the more recent QMC QCapri- 



otti et al., 1999), series-expansion (Zheng et al., 2006c[), 



and density-matrix renormalization group studies IWhite 



|and Chernyshev||2007[ ) confirmed the stability of the 120 c 



spin-structure for the spin-^ case and yielded (S) 



0.20 



for the ordered moment. 

The series-expansion results for the excitation spec- 
trum of the spin-^ triangular-lattice antiferromagnet 



(Zheng et al., 2006a) deviate substantially from the har- 
monic spin-wave theory with the overall band-narrowing 
by ~ 50% and additional rotonlike minima. These effects 
were shown to result from strong magnon-magnon inter- 



actions by |Starykh et al.|fl2006| ) and |Zheng etal] ( |2006c[ ) , 
who demonstrated that the first-order T/S renormaliza- 
tion gives a spectrum in qualitative agreement with the 
series-expansion data. The importance of decays in the 
magnon spectrum of the triangular-lattice antiferromag- 
net w as pointed out by | Chernyshev and Zhitomirsky | 
([2006). 



We now proceed with a brief description of the spin- 
wave results for the dynamics of the triangular-lattice 
antiferromagnet, refer r ing int rested readers to |Cherny-| 
|shev and Zhitomirsky| ( |2009a[ ) for a detailed review. We 
consider the Heisenberg Hamiltonian (38) with the ex- 
change interaction between the nearest-neighbor sites 
J(r) = J5| r n. In this case, the Fourier transform of 
the exchange matrix is J q = 6J7 q with 7 q = ^[cosq x + 
2 cos(q x /2) cos(y/3q y /2)]. The ordering wave vector of 
the 120° spin structure Q = (47r/3,0) is commensurate 
with the reciprocal lattice and is located at the K point 
in the corner of the Brillouin zone; see Fig. [13] 

The magnon dispersion in the harmonic approximation 

is 



£k = ZJS^d - 7 k)(l + 2 7k ) 



(42) 



In full agreement with the general picture outlined above, 
the excitation spectrum has three Goldstone modes at 
k = 0, ±Q with distinct spin- wave velocities cq and cq. 
In the harmonic approximation their ratio is cq/cq = 
y/2, with the leading 1/5 renormalizatio n reducing it 
to ^ Cp /cq 
al. 



1994) 



1.156 in the S = 1/2 case (Chubukov et 



Another notable feature of the dispersion in 
Eq. (|42|) is the varying convexity of the acoustic branch 
in the vicinity of the Q points 



£Q+k ~ c Q k(l - a^k) , 



cos3(^ , 



(43) 



which changes its sign depending on the direction of k. 

Using the harmonic expression for we obtain the 
magnon decay region of the triangular-lattice antiferro- 
magnet as given by the shaded hexagram in Fig. 13 













'f r 







FIG. 13 (color online). Brillouin zone of the triangular lat- 
tice. The magnon decay region is shown as a shaded area. 
Solid and dashed lines indicate saddle-point van Hove singu- 
larities of the two-magnon continuum. 
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The decay threshold boundary is determined solely by 
the emission of an acoustic magnon with q — ^ ±Q; see 
Eq. (17), shown in Fig. 13 by dotted lines. As anticipated 



to magnon interactions is even stronger and exceeds 50% 



above, the decay region includes an extended neighbor- 
hood of the T point. Interestingly, the decay region also 
extends all the way to the K, K' points (= =LQ). This is 
a consequence of varying convexity of the corresponding 
acoustic branch (43) together with the commensurabil- 



ity of the ordering wave vector 3Q = G. As a result, 
the kinematic conditions are satisfied for a decay from 
the steeper side of the energy cone with k — >• Q into two 
magnons with q, q' — >> — Q. Thus, magnons near K or 
K' points are unstable only in a certain range of angles, 
in agreement with Fig. [l3j 

In clear contrast with the field-induced decays (see 
Fig. [5| , the decay region of the triangular-lattice antifer- 
romagnet is not extended by the three-magnon decays. 
This implies that the cascade multimagnon decays do 
not have a tendency to spread over the Brillouin zone 
and that the decay threshold boundary remains sharp. 
The boundary can change due to renormalization of the 
harmonic dispersion (42), but the corresponding shifts 



appear to be no t very significant ( |Chernyshev and Zhit- 
omirsky [ |^)^). 



Last, the solid and dashed lines in Fig. 13 show the 
threshold boundaries for the emission of two equivalent 
and inequivalent magnons, respectively; see Sec. |III.A| 
Both sets of contours lie completely within the decay 
region and thus correspond to the saddle points of the 
two-magnon continuum, which result in the logarithmic 



peaks in the magnon decay rate Tk; see Sec. III.C 



The 1/S renormalization of the magnon spectrum of 
the triangular-lattice antiferromagnet is calculated in the 
same way as in the case of the square-lattice antiferro- 
magnet, discussed in Sec. IV The minor technical sim- 
plification is in the absence of the angle renormaliza- 
tion in the present case due to the high symmetry of 
the 120° state. Aside from this detail, the renormalized 
"on-shell" spectrum is given by Eq. ( [35] ) with the self- 
energies E a (k, £k) an d £&(k, £k) [ see Figs. [TJa,b)], and 
the Hartree-Fock correction from the quartic terms (30). 



All relevant expressions can b e found in|Chernyshev and 



Zhitomirskyl ( |2009a| ); see also Starykh et al. (2006) for a 
somewhat different procedure of deriving the 1/S renor- 
malization. 

The magnon energy and the decay rate obtained in 
the 1/S approximation (35) are shown by dashed lines 
in Fig. [14] along the representative k directions. In ad- 
dition to the singularities, similar to the square-lattice 
case in Fig. [9j there is a significant reduction of the 
magnon bandwidth. This substantial downward energy 
renormalization is due to strong repulsion of one- and 
two-magnon states that are coupled by the cubic anhar- 
monicities. One can expect such interactions to produce 
similarly strong downward shifts of the excitation ener- 
gies in the other spiral antiferromagnets. For example, 
in the anisotropic XY limit of the 5 = 1/2 triangular- 
lattice antiferromagnet, reduction of the bandwidth due 



of the harmonic value (Chernyshev and Zhitomirsky 



2009a), a magnitude unheard of in the collinear antifer- 



romagnets. 

Finally, the magnon damping, the hallmark of magnon 
decays in spiral antiferromagnets, is substantial through- 
out much of the decay region of Fig. [l3j As noted, 
the channel that determines the decay threshold bound- 
ary in the present case is the emission of the acoustic 
magnon. Because of that, the threshold behavior of the 
decay rate near the boundary in Fig. [141 is not sing ular, 
I\ ~ (k — kb) 2 (Chernyshev and Zhitomirsky 2006 ), in 
contrast with the square-lattice case (see Fig. [9|, where 
the threshold singularities are due to minima in the two- 
magnon density of states. 

There is another important difference in the dynamics 
of two-particle decays in these two cases. Magnons emit- 
ted in an elementary decay process in the square- lattice 
antiferromagnet are also unstable. Allowing for a finite 
lifetime of the decay products using self-consistent dress- 
ing of Fig. [TJa) gave us divergence-free magnon damping 
in Sec. |IV| For the triangular-lattice antiferromagnet this 
approach may not work as the emitted magnons are of- 
ten stable. In particular, this is true for the vicinity of 
the logarithmic singularities indicated by solid lines in 
Fig. [13] the saddle-point momenta of the two-magnon 
continuum, which are responsible for the singular be- 
havior, can be shown to fall outside the decay region. 
Hence, at the saddle points, the logarithmic divergence 
of the one-loop diagrams will persist even for the renor- 
malized spectrum and the singularities seem to remain 
essential. Thus, the construction of a renormalized the- 
ory seems to require a summation of an infinite series 
of diagrams that contain leading-order divergences, as 
in the Pit aevskii treatment of the spectrum termination 
problem (Pitaevskii, 1959). Such a procedure is tech- 



nically problematic as it goes beyond the standard 1/S 
spin-wave expansion and allows only for qualitative re- 



sults (Chernyshev and Zhitomirsky 



2009a) 



This seemingly hopeless situation is resolved if we note 
that the logarithmic singularity occurs for a magnon in- 
side the decay region, which means that there is always 
a "background" of the nonsingular decays in Tk- In that 
case, one has to allow for a finite lifetime of the initial 
magnon while magnons created during the decay process 
may remain stable. Technically, the corresponding pro- 
cedure amounts to solving the "off-shell" Dyson equa- 
tion for complex s = — iT^, preserving causality by 
substituting e — >• £* in the argument of the self-energy 
(Chernyshev and Zhitomirsky 2009a). Then the Dyson 



equation is reduced to the system 
e k = e k + Re[E(M% T k = 



-Im[£(k, 



(44) 



The off-shell approach avoids complications related to 
higher-order multiloop diagrams and naturally regular- 
izes the decay singularities in Tk and e^. Results ob- 



tained by numerically solving Eq. (44) are shown in 



Fig. 14 by solid lines. The decay rates remain signifi- 
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FIG. 14 (color online), 
units of J for the spin- 



Magnon energy and decay rate in 
triangular- lattice antiferromagnet. 
Dotted and dashed lines are the harmonic and the 1/S renor- 
malized results, resp ectively. Solid lines are solutions of the 
Dyson equation p4| in the complex plane. 



cant even in the absence of logarithmic peaks reaching 
r k /e k ~0.15. 



VI. DAMPING OF LOW-ENERGY MAGNONS 

The long-wavelength limit is where one usually hopes 
to find universal behavior of the dynamical properties. 
It is also the regime where decay rates can be calcu- 
lated perturbatively because of the smallness of interac- 
tion among low-energy excitations and due to reduction 
of the phase- space volume available for decay processes. 
In this section we discuss a few general scenarios for de- 
cays of the long-wavelength excitations with the goal of 
identifying universal dependencies of r k on k and com- 
plement it with particular examples from the triangular- 
and square-lattice antiferomagnets considered above. 

As discussed in Sec. III.B[ there are several typical 



situations for the overlap of the single-particle acoustic 
branch with the two-particle continuum that allows for 
energy conservation to be fulfilled in the decay process. 
Here we do not repeat this classification and assume that 
such conditions are always satisfied. Generally, the decay 
rate is given by 



r k 



q 



Vk,q| 2 £(£k — £q — £k-q) j 



(45) 



where Vk, q is the decay vertex. Because of energy con- 
servation, the upper limit for the momentum of created 
quasiparticles should be of the order of k. Then, for the 
linear spectrum £ k ~ fc, a naive power counting suggests 



(46) 



where k D comes from the D-dimensional phase space, 
1/k is due to the reduction of that space to the decay 



surface from energy conservation, and is the typi- 
cal amplitude of the cubic vertex on the decay surface. 
We assume that the cubic vertex at small momenta fol- 



lows the standa rd "hydrodynamic" form (Lifshitz and 



Pitaevskii 



[1980J): Vk, q oc \/kqq' , where q r = |k — q|, and 



that in a typical decay process all the momenta are of the 
same order g, q' ~ k. This makes |Vfc| 2 oc k 3 and yields 
a seemingly universal power law for the decay rate of an 
arbitrary acoustic mode: 



r k ~ k D+2 . 



(47) 



For D = 3, T k ~ /c 5 , which matches the well-known re- 
sult fo r the decay of a phonon with a convex dispersion 
curve QBeliaev , 1958), but, as we see shortly, only coinci- 



dent ally. 

In reality, the situation is more delicate and possible 
power-law exponents for T k vary depending on physical 
details. In the case of a single weakly nonlinear acoustic 
branch, relevant to phonons in 4 He, £ k = ck+ak 3 , decays 
are allowed only for the positive curvature a > 0; see 
Sec. III.B A long- wavelength phonon decays into a pair 



of quasiparticles with their momenta directed in a narrow 
solid angle along the direction of the initial momentum. 
The apex angle of the decay cone scales as ~ k such 
that the phase space is k 2D ~ x instead of the naive k D . 
Then the restriction from the energy conservation gives 
1/kO 2 ~ 1/k 3 instead of 1/k as considered previously. 
Altogether, for the case of the cubic upward c urvature of 
the spect rum, the answer is universal [see also Kreisel et 
"^T] ([2008])]: 



L2D-1 



(48) 



k b in 3D ( |Beliaev|[l958 



which indeed yields T k 

and Pitaevskii , |1980|), b ut the answer for 2D is T k 
different from Eq. (47). The 2D result applies to the 



Lifshitz 

f: 



square-lattice antiferromagnet in the strong field consid- 
ered in Sec. HVl 

If several acoustic modes with different velocities are 
present, the fast quasipart icle ca n always decay into two 
slow excitations (see Sec. III.B). This situation is sim- 
pler than in the previous case since the nonlinearity of 
the spectra plays no role. The phase-space factor be- 
comes k D ~ x now, in agreement with the naive considera- 
tion above. Physical realizations of this scenario include 
decays of the longitudina l pho non into two transverse 
ones in solids (Kosevich, 2005) as well as the decay of 
the k —> into two k — >• ±Q magnons in the triangular- 
lattice Heisenberg antiferromagnet considered in Sec. |V| 
a case also pertinent to other spiral antiferromagnets. 
Interaction between phonons in crystals obeys the con- 
ventional asymptotic behavior Vk, q oc \Jkqq r and, conse- 
quently, \Vk\ 2 oc k 3 . Therefore, our initial naive power- 



counting result in Eq. (47) is, actually, valid for them. 
However, in the case of the triangular-lattice antiferro- 
magnet, the result is yet different from that in Eq. (47) 
because the three-magnon vertex is ano malous and scales 
as Vic,Q+q oc (q f — q)\/k/qq l for small k ( [Chernyshev and 
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Zhitomirsky 2009a). For a typical decay process the ini- 
tial and final momenta are of the same order, q' ~ k, 
giving |Vk,Q+q| 2 ~ k instead of k 3 . Altogether, for non- 
collinear Heisenberg antiferromagnets in D dimensions 
this yields 



(49) 



Direct analytic calculation for long-wavelength magnons 
in the triangular antiferromagnet gives a decay rate with 
the correct exponent I\ « 0.79k 2 . 

The decay vertex for the k —> Q magnon has the con- 
ventional scaling VQ+k,-Q+q \fkqq' so the decay prob- 
ability is \Vk\ 2 ~k 3 . However, due to a constraint on the 
angle between k and q, the decay surface in q space is 
a cigar-shaped ellipse with length ~k and width ~k 3 / 2 . 
That causes the restricted phase volume of decays to scale 
as /c( 3jD_5 )/ 2 and leads to a nontrivial k dependence of 
the decay rate: 



fc (3D+l)/2 



(50) 



Numerically, in the triangular antiferromagnet IY 



1.2J&; 7 / 2 along the KT line, in agreement with Eq. (50). 



Away from this direction, damping exhibits an anoma- 
lous angular dependence I\ ~ /c 7//2 /(cos 3(f) 3 / 2 . It di- 
verges upon approaching the decay threshold boundaries 
at (p = ±7r/6 in Fig. 13 and thus requires a self-consistent 
regularization. 

Another case of the singular behavior of the decay 
rate in the long-wavelength limit occurs in the consider- 
ation of the square-lattice antiferromagnet in high fields, 



Sec. [IV] Her e we outline an approa ch to regularizing such 
divergences (|Mourigal et al. 2010). Using the asymptotic 



form of the magnon dispersion in (18), the decay vertex 



V3 ^ (q; k) oc yjkqq' , and the Born expression for the de- 
cay rate in ([7]) we obtain 



r k 



— t^ 2 e x f^k 3 

16tt V 6a 



(51) 



in agreement with Eq. (48) with expressions for c and a 



provided in Eq. (32). An important observation is that 



the nonlinearity of the spectrum a vanishes at the decay 
threshold field H — >> if* and, as a result, the coefficient 
in front of k 3 in Eq. (51) diverges. Such nonanalytic be- 
havior is nothing but the long-wavelength version of the 



decay threshold singularities discussed in Sec. III.C The 
self-consistent regularization ( 37 ) should remain applica- 
ble in this limit. In order to obtain the renormalized k 
dependence of the damping, we substitute the general 
power-law ansatz for the decay rate I\ = /3/c n , and as- 
sume that the damping is much larger than the vanishing 
nonlinearity but still much smaller than the magnon en- 
ergy ak 3 <C j3k n <C ck. In this case, the apex angle of the 
decay cone scales as (p 2 ~ /c n_1 with q ~ fc, which makes 
the decay vertex F 3 (1) (q;k) oc k n ~ 3/2 instead of k 312 . 



The power counting on both sides of Eq. (37) yields a 
unique solution n = 3. Therefore, the Born power law 



Tk = fik 3 is not changed by the self-consistent proce- 
dure. The damping coefficient f3 no longer diverges near 
the decay threshold field because a drops out from the 
self-consistent equation for IY- Instead, j3 exhibits a step- 
like behavior in H — H* , which is in accord with the 2D 
character of the van Hove singularity at the border of the 
two-particle continuum. 

Another interesting asymptotic regime concerns the 
close vicinity of the saturation field with an almost 
parabolic magnon dispersion. As H — » H S1 the velocity 
of the acoustic mode decreases and the condition of weak 
nonlinearity appli es to a pro gressively narrower range of 
momenta k <C 4^1 — H/H s , reducing the range of valid- 
ity of the asymptotic expression (51 ). Outside that k do- 



main one can use the parabolic form for magnon disper- 
sion, £k ~ JSk 2 , to derive another useful asymptotic ex- 
pression for the decay rate. In this regime the decay ver- 
tex is a constant Vg^q; k) ex y/l — H/H s , and the angle 
between the emitted and initial magnons can be large, 
removing restrictions from the decay phase space. Com- 
bining the power-counting arguments following Eq. (46) 
and performing analytical integration in Eq. ^ we ob- 
tain the universal momentum-independent expression 



r k w 16 1 



H 



(52) 



which is valid for 4^/1 - H/H s < k < 1. 



VII. MAGNON DECAY IN QUANTUM SPIN LIQUIDS 

Magnetic materials with quantum-disordered or spin- 
liquid-like ground states have attracted significant atten- 
tion in the past two de cades. The original idea goes 



back to |Anderson| ([1973]) , who suggested the resonating 
valence-bond singlet state as a candidate for a quantum- 
disordered ground state of the spin- \ triangular antiferro- 
magnet. Since then, various types of spin-liquid ground 
states have been discussed and classified; see |Dagotto 



and Rice (1996), Balents (2010), and Lhuillier and Mis- 



guich (2011). On the experimental side, well-known ex- 



amples inclu de the spin-Peierls compound C uGe03 ( Hase 
, the dimer systems CssC^BrQ ( Leuenberger 



et al. 



et al. 



1993 



spin antiferromagnetic chains (Regnault et al. t 1994), and 



1983] ) and TlCuCl 3 ( |Takatsu et a/l|1997[), int eger- 



many others. A common property of all these materials 
is the gap in the excitation spectrum, which separates 
a spin-singlet ground state from the triplet of magnetic 
excitations with 5 = 1. For that reason they are also 
referred to as quantum spin-gap magnets. The spin- spin 
correlations remain short ranged down to zero tempera- 
ture with the correlation length inversely proportional to 
the gap £ ~ 1/A. Because of the three fold degeneracy, 
magnons in quantum spin-gap materials are often called 
triplons. 

The spin-dimer magnets offer a particularly wide va- 
riety of quantum-disordered states regarding the energy 
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scale of the spin gap and different effective dimensionali- 
ties. They became model systems for studying quantum 



phase transitions (Ruegg et a/., 2008 Sachdev 2008) and 



the field- induced Bose-E i nstein condensation of ma gnons 
dGiamarchi and Tsvelik| |1999| |Nikuni et a/.[ [2000| ). 



For 

an overview of this subflel d of quantum magneti s m, se e 
Grenier and Ziman (20071) and iGiamarchi et al. (2008). 



Here we focus on the role of spontaneous magnon decays 
in the dynamics of quantum spin-gap magnets. Thereby, 
we exclude gapless ID spin-liquid states with fractional- 
ized spinon excitations as well as more exotic forms of 
spin liquids in higher dimensions t hat have yet to find a 
definitive experimental realization ( jBalents , 2010). 

The theoretical attention paid to magnon decays 
is motivated by the inelastic neutron experiments in 
two organic quantum spin-gap mate rials p iperazinium 
hexachlordicuprate (PHCC) ( [Stone et al\ |2006|> and 
(CH3)2CHNH 3 CuCl 3 _ called IPA-CuCl 3 QMasuda et al. 



2006 



Zheludev et al. 



2007) in which sharp transforma- 



tion of triplet excitations upon crossing into the two- 
magnon continuum was observed. In PHCC, the magnon 
lifetime is seen rapidly decreasing once inside the con- 
tinuum, whereas in IPA-CUCI3 not even a trace of the 
single-particle peak is detected beyond the crossing. Such 
an end point was dubbed the termination point of the 
magnon branch by analogy with a similar ef fect in the 



superfluid He (Lifshitz and Pitaevskii 1980). The ori- 



gin of this drastic effect was linked to the nonanalytic 
behavior of the magnon self-e nergy due to two-particle 
decays by Zhitomirsky ([2006 ). The predicted behavior 



was confirmed by numerical simulations (jFischer et al. 



2010 Zheng et a l. 2006b) and by an analytical model 
calculation ( |Bibikov| |2007| ). In this section we give an 
overview of theoretical results on two-magnon decays 
in quantum-gap antiferromagnets. We discuss a mi- 
croscopic approach to spin-dimer magnets based on the 
bosonic representation of spins in terms of bond oper- 
ators in Sec. |VII.A| This is followed by an analysis of 
the termination point in the magnon dispersion curve; 
see Sec. |VII.B| Discussion of the related problem of the 
decay of longitudinal magnons in weakly ordered antifer- 
romagnets is included in Sec. |VII.C| 



A. Bond-operator formalism 
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FIG. 15 (color online). Geometry of exchange bonds for (a) 
symmetric spin ladder and (b) asymmetric ladder with extra 
diagonal coupling. 



the rungs of the ladder. In this limit, properties of the 
spin ladder can be calculated using the bond-operator 
representation of spins proposed independently in sev- 



eral works ( Chubukov 
iBhattl 11990 ~ 



1989 Harris, 1973 Sachdev and 



Two spins S = 1/2 on the same rung can form one 
singlet I s) and three triplet states: 



\s) = 
l*o> 



1 

71 
1 

: 71 



(in)-ut», |i+> = in>, 
(in) + lit)), i*-) = iw>. (54) 



In zero field it is more convenient to use the vector basis 
\t x ) = - \t+))/y/2, \t y ) = i ( \ t-) + |t+»A/2, and 

\tz) = \to)'i see Sachdev and Bhatt ( |1990[ ). The represen- 



tation of the original spin operators in the dimer basis 
is 



°1,2 



±7 ) (s\ + tls)--e^tlt 1 , (55) 



where we introduced singlet and triplet tf a operators 
that create four physical states out of a fictitious vacuum 
|0). Bosonic commutation relations for s and t a together 
with the constraint s^s + t\ t t a = 1 ensure the standard 
SU(2) algebra for spins. 



In the bond-operator representation (55), the spin 
Hamiltonian (53) takes the form T~L = Ho + ^2 + 7^4 > 



One of the most studied models of quantum spin-dimer where H and U 2 are quadratic in t a : 



systems is the spin- 
Hamiltonian is 



ladder shown in Fig. 15 'a). Its 



H — J± ^ Sii • + J ^ S ni • S ni+ i , (53) 

i i,n=l,2 

where J± is an antiferromagnetic exchange along rungs 
and J is a coupling along the legs of the ladder. Such 
a model spin Hamiltonian is realized in a few magnetic 
compounds ; see, f or example, Notbohm et al. ( 2007[ ) and 
Schmidiger (2011). For J± > J, the ground state can be 
approximated by a product of singlet dimers occupying 



Ho 



t ] t 



)■ 



2 



t 



i+1 m+1 ai 



' S \ S \+l^ai+l^ai 



(56) 
h.c.) , 



whereas H4 includes four-triplon terms. Next one needs 
to enforce the local constraint on the total number of sin- 
glets and triplets on each dimer. This can be achieved 
using various approximate schemes. One of them con- 
sists of condensing singlets Si,s| — >• (5), and fixing the 
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boson density with the help of the chemical potential 
flSachdev and Bhattj [1990 ) . Alternatively, the constraint 
is resolved in Holstein-PrimakofT-like fashion, s = = 

— tat a , and the square roots are expanded in pow- 
ers of the triplon density (Chubukov, 1989). Although 



the final form of the effective Hamiltonian is somewhat 
different in the two approaches, the lowest-order magnon- 
magnon interaction contains only four-particle terms. 

Suppose now that the crystal space group allows ad- 
ditional asymmetric diagonal coupling Jd between rungs 
of the ladder [see Figfl5^b)]: 



2i+l • 



(57) 



The bond-operator representation for Hd yields a three- 
magnon interaction term 



= y ^)E ea/37 ( sin ^4 1 XAs7+ h ' c ')' (58) 



t ,t 



which has the structure anticipated in Sec. II. D 

The difference in the magnon interaction between the 
two models [see Figs. 15 ^a) and (b)] has a symmetry 
origin. Three-magnon interactions are prohibited in a 
symmetric spin ladder because singlet states are anti- 
symmetric under permutation of two spins while triplet 
states have even parity. However, if the reflection sym- 
metry between the ladder legs is broken, coupling be- 
tween one- and two-triplet sectors (58) is also allowed. 



Three-magnon interactions exist in many spin-dimer sys- 
tems including the bond-alternating spin chain and var- 



ious 2D and 3D arrays of dimers (Matsumoto et al. 



2004 Sachdev and Bhatt 



1990). Cubic terms have also 



been discussed for the Haldane-chain antiferromagnets 
flKolezhuk and Sachdev[|2006||Zaliznyak et a/.||20Qll ). Fi- 



nally, the geometry of exchange bonds in the spin-dimer 
magnet IPA-CUCI3, which will be discussed in the next 



section, has precisely the structure shown in Fig. 15 'b) 
(Fischer et all 120111 iMasuda et all 12006). 



B. Termination point of the magnon branch 

We now investigate the effects of three-magnon inter- 
actions for the generic energy spectrum of a quantum 
spin-gap magnet illustrated in Fig. 16 An important 
feature of the magnon dispersion is a minimum at a fi- 
nite momentum due to antiferromagnetic correlations in 
the system. For simplicity, we consider only the most 
common case of a single minimum at a commensurate 
wave vector Q. The overlap between the single-magnon 
branch and two-particle continuum occurs if the magnon 
bandwidth W = O(J) exceeds the gap A = (D(J±); see 



Fig. [16} Then, naturally, the decay region is the range of 
momenta centered around k = and magnons that are 
created in a decay process are farther down the dispersion 
curve toward Q and — Q. 
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FIG. 16 (color online). Structure of the energy spectrum of 
a quantum spin-gap magnet. The thin solid line is the bare 
triplon energy £k, the full solid line is the renormalized disper- 
sion £k, and the shaded region is the two-particle continuum. 
The decay threshold is marked by k c and e c . 



Following the kinematic analysis of Sec. |III.A[ in the 
absence of gapless Goldstone modes, the only two chan- 
nels that may define the boundary of the decay region are 
the decays into magnons with equal momenta or equal ve- 
locities [cases (i) and (iv), respectively]. Considering the 
typical situation of a modest overlap with the continuum, 
i.e., when the magnon energy e c > 2 A, one can show that 
the prevailing channel is the former channel of equal mo- 
menta £k c = 2£(k c +G)/2? where G is a reciprocal lattice 
vector. The threshold boundary for the field-in duce d de- 
cays in square-lattice antiferromagnet (see Sec. IV), is of 
the same type, yet the kinematics is quite different here 
as the emitted magnons are outside of the decay region 
and no cascadelike smearing of the boundary is possi- 
ble. Thus, we encounter a kinematic situation distinct 
from the cases treated in Sees. [TV| and \V\ which requires 
detailed consideration of the spectrum near the decay 
boundary. As shown next, the presence of nonanalytic 
self-energy leads to dramatic qualitative changes of the 
magnon spectrum near the threshold. 

According to the preceding discussion, the three- 
magnon interaction has the form 

^3 = ^^(q^-^kje^^^^+h.c.). (59) 

k,q 

Although due to conservation of the total spin the 
magnon with polarization 7 can decay only into magnons 
with different polarizations, as enforced by the antisym- 
metric tensor e a/37 , this has no implications for the decays 
as long as the three branches remain degenerate. More 
important is the antisymmetric structure of the magnon 
vertex under the permutation of polarizations a and f3 
and, consequently, under the interchange of momenta of 
the decay products ^(q, q';k) = — T^(q', q; k). Since 
at the decay boundary the emitted magnons have equal 
momenta, the three-magnon vertex (59) vanishes at the 
threshold Va(q*,q*;k c ) = 0. Then the main difference 



from the analysis of the decay singularities in Sec. III.C 
is the effective increase of dimensionality due to further 



22 



expansion of the decay vertex in q = q 
the vicinity of the boundary: 



S(k,w) 



"Si 



JO 



q D+1 dq 



Ag 2 +v 2 -Ak- Auj-iO 



q* needed in 



(60) 



where D is the dimension, Auj = uj — e c , Ak = k — 
k c , V2 is the velocity of emitted magnons, A = O (J) 
is a constant, and g% = O(Jd) is an effective coupling 
constant. Because of the vanishing vertex, the self-energy 
(60) is not singular at the decay boundary Acj, |Ak| — » 
even in D = 1, but it is still nonanalytic. The higher- 
order corrections can be shown to leave unchanged the 
functional form of the nonanalyticity obtained in the one- 
loop approximation (60). Next we absorb the regular 



part of E(k, uj) into the magnon energy while the 
renormalized spectrum will be obtained from the Dyson 
equation with the nonanalytic part E(k, uj). We focus 
on the ID case, mentioning a few results for the higher- 
dimensional systems in the end. 

We begin with the region k > k c where spontaneous 
magnon decays are forbidden; see Fig. fT6] In this case, 



integration in Eq. (60) yields the following nonanalytic 
self-energy: 



£(k,u) Xy / v 2 Ak - Auj , Auj < v 2 Ak (61) 

with A = 7n?!/2A 3 / 2 = 0(J 2 /J 3 / 2 ). Close to the decay 
boundary the nonanalytic E(fc,u;) dominates the regular 
part of G~ 1 (k^uj) ~ Ak,Auj and completely determines 
the behavior of the quasiparticle pole. Solving the Dyson 



equation with Eq. (6l[ ) yields the renorm alized spectrum 
sketched in Fig. 16 (Zhitomirsky| 2006). Far away from 



the decay boundary the renormalized magnon energy Sk 
follows bare dispersion Sk with the negative velocity v\ = 
dek/dk < 0. On the other hand, for k —> k c the renormal- 
ized slope must change its sign to det/dk = v 2 > as the 
minimum of the two-particle continuum at k c is formed 
by magnons from the positive slope of This change of 
the magnon velocity takes place in the crossover region 



Ak ~ X 2 /(v 2 - vx) 



0(J^/J 2 ) and the overall shape of 



€k resembles the avoided crossing between the magnon 
branch e k and the continuum boundary. However, due 
to the nonanalytic nature of the self-energy the behavior 
of the spectrum is markedly different from the avoided 
crossing: the quasiparticle weight Z k changes drastically 
in the crossover regime and vanishes at the intersection 
of the magnon branch with the two-magnon continuum. 
Such a suppression of the quasiparticle peak is a hall- 
mark of the termination point in the energy spectrum, 
very similar to that in 4 He ( |Pitaevskii| |1959[ ). In a dif- 
ferent context, similar distortion of the excitation curve 
Sk in the vi cinity of the continuum boundary was also 
observed by |Coldea et al. (2010). 

For k < fc c , inside the continuum where spontaneous 



decays are allowed, the nonanalytic self-energy (60) is 
now purely imaginary: 



E(fc, uj) = -i\y/\v 2 Ak - Auj\ , Auj > v 2 Ak . (62) 
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FIG. 17 (color online). Time-of-flight inelastic neutron- 
scattering d ata measured in IPA-CuCl3 at T — 1.5 K ( |Masuda| 
et al\ 2006 ) . The arrow indicates the termination point of the 



magnon branch. From Andrey Zheludev. 



Within the on-shell approximation, uj 
the following magnon decay rate: 



this leads to 



T k ^ A v / (^2 - vi)\k - k c \ . 



(63) 



The s ame result was obtained by[Kolezhuk and Sachdev| 
(2006). A more rigorous off-shell solution of the Dyson 



equation in the complex plane changes the perturbative 
result of Eq. ^ in the Afc vicinity [~ C(Jj/J 2 
the decay boundary to T k oc \k — k c \ 2 ^ 3 . However 



of 



more 

importantly, the quasiparticle weight remains suppressed 
in this region, in the same way as on the other side of the 
termination point. The quasiparticle peak may reappear 
farther away from the decay boundary depending on the 
strength of the three-magnon interaction. 

We also note the common misconception regarding the 
absence of a pole beyon d the termination p oint inside the 
two-particle continuum ( |Stone et a/^|2006| ). According to 
our analysis, such a pole is always present on either side 
of the termination point, but its quasiparticle residue is 
vanishingly small in the vicinity of that point. 

The effect of spontaneous decays in the energy spectra 



Zheng et al ( 


2006b) 


Fischer et al. 





bond-alternating spin-^ chain by 
and for the asymmetric ladder by 
For both models, the single-particle branch merges with 
the lower edge o f the continuum, gr adually losing its 
spectral weight. |Fischer et al.\ ([2010 ) also found some 
indications of a reemerging quasiparticle peak inside the 
continuum for small values of the three-magnon coupling 

J A. 

A remarkable illustration of the theory is provided 
by the energy spec tra measured in the s pin-ladder com- 
pound IPA-CuCl 3 flMasuda efal] |2006| ). This material 



consists of 2D arrays of weakly coupled asymmetric S ■ 
1/2 Heisenberg ladders [see Fig. [l5]fo)], with the coupling 
constants J± = 2.9 meV, J = 1.2 meV, Jd = —2.3 meV, 
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and a small interladder exchange J id 
cher et al\ |2Q11 ). The resolution-limited peaks cor- 



0.3 meV (Fis- threshold singularities essential. 



responding to the magnon band with a minimum at 
Q = 0.5 (units of 2tt /a) and the energy gap A = 1.2 meV 
are observed in IPA-CUCI3; see Fig. [17 



magnon branch a bruptly disapp e ars at k c 
their data for 



However, the 
« 0.8. Using 



Masuda et al. (2006) determined that 



k c indeed corresponds to the intersection point of with 
the continuum. Thus, the single-magnon branch in IPA- 
CUCI3 terminates at the boundary of the two-magnon 
continuum, as in the superfluid 4 He, where the quasi- 
particle branch te rminates at the threshold of two-roton 
decays flGfyde et al\ |1998| [Pitaevskii| |1959[ ). 

A gradual broadening of the quasiparticle peak by two- 
magnon decays without its complete disappearan ce was 
observed in the 2D dimer organic material PHCC (Stone 
et al.\ 2006]) and in the q uasi- ID antiferromagnet CSN1CI3 
flZaliznyak et a/.l[200l|. The difference from IPA-CuCl 3 



discussed above is mainly due to the higher dimensional- 
ity of these two materials. In the case of PHCC, the ex- 
perimental data are reasonably well described by Fermi's 
golden-rule formula 



A|k-k c |^ 2 , 



(64) 



with D = 2 (Kolezhuk and Sachdev, 2006). Beyond 



the on-shell approximation, the triplon self-energy (60) 



in D > 1 gives much weaker nonanalyticity. For exam- 
ple, in 2D, E(k,cj) ~ Auj ln(A/Acj), and the suppression 
of the quasiparticle peak in PHCC should occur only in 
an exponentially small region near the crossing point, 
justifying the validity of the perturbative result (64). 



Thus far, we discussed spontaneous magnon decays in 
quantum spin-gap systems under the "standard" condi- 
tion of the cubic vertex vanishing at the decay threshold 
boundary. The decay-related singularities may be sub- 
stantially enhanced, however, if the vertex is finite. This 
takes place, for example, when the threshold channel cor- 
responds to decays into magnons with different momenta. 
Another mechanism is spin anisotropy, which violates 
spin conservation and removes the odd parity of the cu- 
bic vertex ( 59 ) with respect to spin indices and momenta, 



making the decay vertex finite at the continuum bound- 
ary and also lifting degeneracy of the triplon bands. This 
mechanism might be relevant to the spin-1 alternating 
chain antiferromagnet Ni(CgH24N4)(N02)C104 (called 
NTENP), whose unusual spin dynamics in external mag- 
netic field may signify s pontaneous two-magno n decays 
(Hagiwara et al. 2005 Regnault et al, 2006). Finite 



magnetic fields also produce intriguing changes in the 
dynamical properties of IPA-CUCI3 (|Garlea et al. , 2007 



Zheludev et al. 2007). Theoretical explanation of these 



experimental findings remains an open problem. 

Finally we note that the occurrence of the termination- 
like behavior in the spectra of quantum spin-gap magnets 
is by no means unique. It can also take place in a vari 
ety of spin-anisotropic antiferromagnets ( Chernyshev and 



Zhitomirsky| |2009a[ ) , in which the anisotropy gap stabi- 
lizes the decay products and makes the corresponding 



C. Decay of longitudinal magnons 

Our final example of spontaneous magnon de- 
cays concerns weakly-ordered collinear antiferromag- 
nets, which occupy an intermediate place between 
quantum-disordered and well-ordered semiclassical mag- 
nets. Apart from the conventional spin waves, excita- 
tion spectrum of the weakly-ordered spin system may 
display a longitudinal mode that corresponds to fluctua- 
tions of the magnitude of the order parameter ( |Affleck| 
1989). The presence of this mode can be easily under- 
stood by considering a transition from the spin-liquid 
phase int o an ordered state under an external pressure 
P ( |Matsumoto et al\ |2004[ ) . In the quantum-disordered 
state, up to the critical pressure P c where the gap closes, 
three magnon modes are degenerate. Above the quan- 
tum transition, P > P c , two of these modes transform 
into the Goldstone modes of the ordered magnetic struc- 
ture, while the remaining third excitation acquires a finite 
gap. At this point one may ask how the above scenario 
could at all agree with the conventional spin-wave pic- 
ture, which accounts only for the two magnon modes. 
The apparent contradiction is resolved if we notice that 
the longitudinal mode is always unstable with respect to 
decay into a pair of trans verse spin waves and thu s has 
an intrinsic linewidth (Affleck and Wellman, 1992). Far- 



ther inside the parameter domain of the ordered phase, 
the longitudinal excitation transforms into two-magnon 
resonance and eventually disappears. 

The pressure-induced quantum critical point is re- 
alized, for instance, in the spin-dimer antiferromagnet 
TICUCI3 ( Oosawa et al. 2003 ) . The low critical pressure 



P c ~ 1 kbar enables the use of inelastic neutron measure- 
ments to explore the phases on both sides of the critical 
point. The existence of the gapped longitudinal mode 
in the antiferromagnetic phase of TICUCI3 was demon- 
strated by Ruegg et al ( 2008 ) , who concluded that such 
a mode is "critically well defined" because of its large 
broadening. Recently, the lifetime of longitudinal excita- 
tions in TICUCI3 was calculated by |Kulik and Sushkov 
(2011). In the following, we describe basic aspects of 



spontaneous magnon decays in weakly ordered collinear 
antiferromagnets using the example of TICUCI3. 

The harmonic theory of magnetic excitations in both 
dis ordered and ordered pha ses of TICUCI3 was developed 
by |Matsumoto et al. (2004) who used the bond-operator 
formalism. At the critical pressure P = P c the gap of 
the triplon branch closes at the commensurate ordering 
wave vector Q. No spontaneous decays of the low-energy 
magnons exist up to this point. At P > P c one member 
of the degenerate triplet condenses, U z —> (t z )e l( ^' ri , and 
the collinear magnetic structure develops along the corre- 
sponding direction. The degeneracy of the triplon bands 
is removed in the ordered phase. Transverse excitations 
t XiV stay degenerate with each other and remain gap- 
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less with the energy £q +k ~ cfc, whereas the fluctuating 
part of the condensed magnon t z = t z — (t z ) represents 
the longitudinal mode with the energy and the gap 



Although the cubic vertex (59) may be absent in the 
disordered state at P < P c , the triplet condensation in 
the ordered phase always produces a three-particle in- 
teraction involving longitudinal fluctuations; see Eq. @: 



6 ak 4 6 ak 3 6 /3k 2 6 /3ki 



(tz)t 



2k 3 6 /3k 2 



t /3ki . Thus, the symme- 



try ensures that the longitudinal magnon can decay into 
a pair of transverse magnons j3 = x,y. Both energy and 
momentum conservation are also trivially satisfied for the 
decay of the gapped mode in the vicinity of the ordering 
wave vector Q. The asymptotic for m of the decay vertex 
i S y 3 (q,q';k) ~ 

|2011| ). In the bond-operator forma 
denominator in V3 with £q +q 



(Kulik and Sushkov 
ism, the anomalous 
cq appears from the co- 
efficients of the Bogolyubov transformation. The decay 
rate of the longitudinal magnon at k = Q calculated from 
Eq. is 



Q 



7Ai 



(65) 



where the final relation follows from the mean-field re- 
sult for the expectation value of the condensate fraction 
(t z ) oc y/P — P c . The theoretical estimate of the damp- 
ing coefficient 7 = 0.2 is not far from the experimental 
value 7 ~ 0.15 and can be further improved upon by tak- 



ing into accoun t spin anisotropy in TICUCI3 (Kulik and 
|Sushkov| |2011[ ). Thus, the decay rate of the longitudi- 
nal mode is only parametrically small compared to its 
energy and whether it remains well-defined depends on 
microscopic parameters of the system. 



VIII. SUMMARY AND OUTLOOK 

We presented an overview of the phenomenon of spon- 
taneous magnon decays in quantum antiferromagnets 
and discussed symmetry and kinematic conditions that 
are necessary for it. While no exhaustive list of magnets 
that favor quasiparticle decays can be given, the offered 
examples outline the principal classes of spin systems, in 
which the effect may occur. Generally, the noncollinear 
spin structures in the systems where the symmetry is bro- 
ken spontaneously from the high-symmetry spin Hamil- 
tonian belong to one of such classes. The gapped quan- 
tum spin liquids is yet another class and a proximity to 
a quantum critical point can be another common situa- 
tion in which magnetic excitations may be intrinsically 
damped. 

The hallmarks of the decay-induced anomalies in 
magnon spectra include substantial broadening of 
magnon peaks in a large part of the Brillouin zone, 
termination points, and non-Lorentzian features in the 
structure factor, as well as strong deviations of the spec- 
trum from the harmonic theory, especially for S = 1/2 



systems. With improvement in the resolution of the 
neutron-scattering experiments expected in the near fu- 
ture, these effects should become readily observable. Be- 
sides their effect on the dynamical properties, magnon 
decays may also play a role in the critical behavior in the 
vicinity of quantum critical points (Fritz et al, 2011). 



We also highlighted the role of singularities, which are 
generally concomitant with decays and are related to the 
crossing of the single-particle branch with the surface 
of van Hove singularities in the two-particle continuum. 
Regularization of such singularities represents a theoret- 
ical challenge that is largely unfamiliar in the traditional 
spin-wave theory and requires the use of self-consistent 
schemes exemplified in this work. While most of the sys- 
tems considered here are ID or 2D, for which the effects 
of spontaneous magnon decays are amplified, the strong 
effects due to spontaneous decays are not restricted to 
low-dimensional systems, as shown in Sec. |IV| 

The direct experimental evidence of magnon decays in 



low-dimensional spin-dimer systems discussed in Sec. |VII| 
is unequivocal and progress is being made toward similar 
observations in the noncollinear and field-induced set- 
tings. For the latter case, several materials with suffi- 
ciently small interactions that allow magnetic fields to fit 
into the window accessible by neutron-scattering experi- 
ments have recently been synthesized and the first indi- 
cation of field-induced magnon decays has been found in 
the spin- 1 squa re-lattice antiferromagnet Ba2MnGe207 



( |Masuda etal] [2010). As noted in Sec. [TV] the BEC 



magnets are also very promising for the observation of 
the field-induced decays in the vicinity of their lower 
quantum-critical point, which is readily accessible in a 
number of them. The experimental search for sponta- 
neous magnon decays in noncollinear magnets is also on- 
going, with several S = 1/2 t riangular-lattice an tifer- 
romagnets, su ch as Cs2CuBr4 (|Tsujii et al. 2007) and 
Ba3CoSb20g (Shirata et al. , |2012 ), and also the large- S 



weakly anisotropic noncollinear antiferromagnets (Ishii et 



al. 2011 Toth et al., 2011) being most promising. 



Besides improvements in conventional neutron scatter- 
ing, there has been a recent dramatic improvement in 
neutron- spin-echo spectroscopy technique; see Pap pas et 
al\ (12006). The energy resolution of this method can 



reac 



1 the /ieV range compared to the meV resolution 
for triple-axis spectroscopy. Successful measurements 
of magnon lifetimes by the neutron-spin-echo technique 
were r eport ed by Bayrak ci et al. ( 2006 ), |Nafradi et al.\ 
(2011), and Chernyshev et al. ( |2012| ) and more of such 
studies are forthcoming. This development may offer an 
opportunity to observe a variety of new physical effects 
including spontaneous magnon decays. 

Finite magnon lifetimes can also be detected by exper- 
imental probes other than the neutron scattering. One 
such example is the decay-related smearing observed by 
electron spin resonance in the bond-alternating chain an- 
tiferromagnet NTENP QGlazkov et~aT\ |2010|). Another 



promising experi mental technique is resonant inelastic 
x-ray scattering (Ament et al., 2011), which is able to 
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detect excitation broadening in large-exchange spin sys- 
tems. The magnetic component of thermal conductiv- 
ity should also be strongly influenced by spontaneous 
magnon decays (Kohama et al t 2011). 
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